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Abstract 

In these lectures I will review the approach to glasses based on the replica formalism. Many 
of the physical ideas are very similar to those of older approaches. The replica approach has the 
advantage of describing in an unified setting both the behaviour near the dynamic transition (mode 
coupling transition) and near the equilibrium transition (Kauzman transition) that is present in 
fragile glasses. The replica method may be used to solve simple mean field models, providing explicit 
examples of systems that may be studied analytically in great details and behave similarly to the 
experiments. Finally, using the replica formalism, it is possible to do analytic explicit computations 
of the properties of realistic models of glasses and the results are in reasonable agreement with 
numerical simulations. 
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1 Introduction 



1.1 General considerations 

In these recent years many progresses in the study of glasses have been done using the replica formalism. 
There are many indications that, if we would follow the evolution of a glass at a microscopical level, 
we would discover that at low temperatures the glass freezes in an equilibrium (or quasi equilibrium) 
configuration that is highly non-unique. This essential non-uniqueness of the ground state is present 
in many others systems where the energy landscape is highly corrugated: e.g. it is widely believed 
to be present in spin glasses, i.e. magnetic systems with a random mixture of ferromagnetic and 
antiferromagnetic interactions [1, 2, 3, 4]. This property is responsible of the peculiar behaviour of 
glassy systems and, at the same time, it is the most difficult to control theoretical. 

The replica approach to spin glasses was developed in the seventies and it was shown to be the 
most effective and sophisticated tool to study the behavior of systems characterized by the presence 
of many equilibrium states. The use of the same replica techniques to study not only spin glasses, 
but also glasses, was not immediate because there was a a strong phycological barrier to be crossed: 
disorder is present in the Hamiltonian for spin glasses, but it not present in the Hamiltonian for glasses 
and it was believed that the presence of disorder was one of the prerequisites for using the replica 
method. However at a certain moment it became clear that disorder in the Hamiltonian was not 
necessary for using the replica techniques and the study of glasses using the replica method started 
to develop rather fast. 

In the picture glasses may freeze in many microscopically different configurations when we decrease 
the temperature. This statement is common to many other approaches [5, 6], however the replica 
approach gives us a panoplia of sophisticated physical and mathematical tools that strongly increase 
our ability to describe, study and compute analytically the properties of glasses. 

These tools have been used to compute analytically in a detailed way the properties of toy models 
for the glass transition. Although theses models are very far from reality (the range of the forces 
is infinite) they display a very rich behaviour [7]: for example there is a soluble mean field model 
without quenched disorder where there is an equilibrium glass-liquid transition (Kauzman transition 
[8]), a dynamical transition[9, 10, 11] (mode coupling transition([12]) and, at an higher temperature, 
a liquid-crystal transition that cannot be seen in the dynamic of the system (starting from the liquid 
phase) unless we cool the system extremely slowly [13]. The existence of these soluble models is very 
precious to us; they provide a testing ground of new physical ideas, concepts, approximation schemes 
that are eventually used in more realistic cases. 

The aim of these lectures is to present an introduction to the replica tools, to describe some of the 
technical details and to stress the physical ideas. The amount of work that has been done in the field 
is extremely large and here I have been forced to concentrate my attention only on a few points. 

1.2 Classiness, metastability and hysteresis 

An essential feature of a glass system at the microscopic level is the existence of a corrugated free 
energy landscape. One may wonder which are the macroscopic counterparts of this property. A very 
important consequence is the presence of metastability in an open region of parameter space, a new 
and unusual phenomenon that can be experimentally studied in a carefully way. 

Let us describe a non-glassy system where we have metastability. The simplest case is a system 
that undergoes a first order phase transition when we change a parameter. When the first order 
transition happens by changing the temperature, if we cool the systems sufficiently slowly, the high 
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Figure 1: An artistic view of the free energy of a system with corrugated free energy landscape as 
function of the configuration space. 



temperature phase survives also below the critical temperature up to the spinodal temperature. 

In order to present a familiar example I will consider a system where the control parameter is the 
magnetic field h: the simplest case is the ferromagnetic Ising model. At low temperature the equilib- 
rium magnetization m{h) is given by m{h) = nig sign{h) + 0{h) for small h {nis being the spontaneous 
magnetization): the magnetization changes discontinuously at /i = in the low temperature phase 
where / 0. 

Let us consider a system that evolves with some kind of local dynamics. If we slowly change the 
magnetic field from positive to negative /i, we enter in a metastable region where the magnetization 
is positive, and the magnetic field is negative. The system remains in this metastable state a quite 
large time, given by t(/i) tx exp(A/|/i|'^), where a = d — 1 [14]. When the observation time is of 
order of t(/i) the system suddenly jumps into the stable state. This phenomenon is quite common: 
generally speaking we always enter into a metastable state when we cross a first order phase transition 
by changing some parameters. 

If we start with the state where m > at /i = and we add a positive magnetic field h at time 0, 
the linear response susceptibility is equal to 

d 

XLR = lim -^Mt, h), (1) 

t^oo an 

m{t, h) being the magnetization at time t. Using general arguments we can show that 



lim > 



{a{t)a{^)Y ^ lim V((a(i)a(0)) - (a(.))(<T(0))). 



(2) 



The linear response susceptibility is not equal to the equilibrium susceptibility that at h exactly 
equal to zero is infinite: 



d_ 

dht 



lim m{t,h) 



h=0 



d_ 

dh 



sign(/i) 



oo. 



h=0 



(3) 
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Indeed 



Xeq{h) = XLR{h) + ms5{h) , 



(4) 



and the two susceptibihties Xeq and xlr differs only at /i = 0. 

This is the usual stuff that is described in books [15]. We claim that in glassy systems the situation 
is different. For example in the case of glassy magnetic systems (e.g. spin glasses) there should be 
an open region in the space of parameters, where, if we change a parameter of the system (e.g. the 
magnetic field h) by an amount A^, we have that xlr 7^ Xeq- We expect that for \h\ < hc{T) we stay 
in the glassy phase ^. 

In this region 



The glassy phase is thus characterized by a non-zero value of Xirr [2]. If we observe the system 
for a time less that r(A/i), the behaviour of the system at a given point of the parameter space 
depends on the previous story of the system and strong hysteresis effects are present. I stress that, 
using the previous definitions, hysteresis and history dependence do not necessary imply glassiness. 
Hysteresis may be present if the time scale for approaching equilibrium is very large (larger than the 
experimental time), but finite. Glassiness implies an equilibration time that is arbitrarily large. In 
other words hysteresis can be explained in terms of finite free energy barriers that may involve only 
a few degrees of freedom; glassiness implies the existence of arbitrarily large barriers that may arise 
only as a collective effect of many degrees of freedom. 

From my point of view the aim of the theoretical study of glasses is to get a theoretical under- 
standing of these effects and to arrive to a qualitative and quantitative control of these systems. 

In these lecture I will address to various aspects of disorder systems using the replica approach. 
The replica formalism is a language for studying different properties of disordered systems: as any 
other language has some advantages and some disadvantages. The replica language has the advantage 
of being very compact and of putting the dirty under the carpet (thais may be also a disadvantage in 
some situations) . In other words you can do computations in a very simple and effective manner, but 
sometimes it hard to understand the meaning of what are you doing. The formalism of replicas and 
replica symmetry breaking seems to be the most adequate to discuss from a thermodynamic point of 
view the situation where Xirr / 0. 

The main physical problem I would like to understand is the characterization of the low temper- 
ature phase in structural glasses and in spin glasses. As you will shall see in these lectures and in 
Cugliandolo's lectures, the main characteristic is the presence of aging for the response[16, 17, 18]. 
Let us consider an aging experiment where the system is cooled at at time 0. 

The response function R{t, tw) is the variation of the magnetization when we add an infinitesimal 
field at time t^. Aging implies that the function i?(f,t^) is not a constant in the region where t and 
tw are both large; indeed we have 



Am{t) = XLRAh for 1<< t « r(A/i), 
Am{t) = Xeq^h for r(A/i) << t, 

where r(A/i) may have a power like behaviour (e.g. r(A/i) cx |A/i|~^). 
It is convenient to define the irreversible susceptibility by 



(5) 



Xeq — Xlr + Xirr ■ 



(6) 



R{t,tyj) = Rs for t « t, 
R{t,tw) = Re for t»t. 



(7) 



^The function hc{T) increases when we decrease the temperature; hc{T) vanishes at the critical point. 
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By definition Rs = xlr and the identification oi Re witli Xeq follows from general arguments that 
will discussed in Cugliandolo's lectures. 

As we shall see later, replica symmetry is broken as soon Rs 7^ Re] we will assume that this 
happens in many glassy systems. This is what is experimentally seen in experiments done by humans 
using a value of ty^ that is much shorter of their life time. These experiments are done for times that 
are much larger (15 or 20 order of magnitude) that the microscopic time and many people do believe 
that aging in the response, i.e. Rg 7^ Re, survive in the limit where the waiting time goes to infinity. 

One could also take the opposite point of view that aging is an artefact coming from doing the 
experiments at too short times and that aging would fades in the limits tw ^ 00. This would be a 
quite different interpretation of the experiments; in these lecture I will stick to the hypothesis that 
aging is present and that Rs remains different from Re also in the infinite time limit. Of course one 
could use experiments to decide which of the two hypothesis is correct, however this must be done 
by extrapolating the experimental data and the discussion would be too long to be presented here: 
moreover, we should remember that quite often, due to practical limitations, the experiments, like the 
God whose siege is at Delphi, neither say nor hide, they hint. 

1.3 The organization of these lectures 

These lectures are organized as follows. 

In section II I will present the simplest model of a glassy system: the random energy model 
(REM) ; I will describe the analytic solution of this model using both a direct approach and the replica 
formalism. 

In the next section I will study models with correlated disorder (e.g. p-spin models); the analytic 
solution of these models is much more complicated than that of the REM and it can be done using 
the replica formalism; it this framework the meaning of spontaneously broken replica symmetry is 
elucidated. 

In section IV I will introduce a key concept of this (and of others) approaches: complexity (aka 
configurational entropy). A particular attention is given to the definition of this quantity: there is 
an intrinsic, albeit small, ambiguity in the value of the complexity that is present in short range 
model, but is absent in infinite range mean field models (if this ambiguity is neglected contradictory 
results may be obtained). Various methods, both analytic and numeric, are introduced to compute 
this fundamental quantity. 

In the next section I will present some general structural properties: equilibrium stochastic stability 
and its dynamic extension that can be used to prove the generalized fluctuation dissipation relations 
that are described in more details in Cugliandolo's lectures. 

In section VI I will present a very short introduction to the physics of structural glasses. 

In the next section I will describe an analytic approach that can be used to do explicit first 
principles computations of the properties of glasses in the low energy phase. Different approximation 
techniques are illustrated. 

In section VIII I will present the results obtained using the techniques of the previous section for 
computing the properties of systems of interacting particles with a realistic potential. These analytic 
results compare in a favorable way to the numerical simulations. 

Finally in the next section one finds a discussion of the open problems. Four appendices dedicated 
to technical problems close these lectures. 
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2 The random energy model 



2.1 The definition of the model 

The random energy model [19, 20] is the simplest model for glassy systems. It has various advantages: 
it is rather simple (its properties may be well understood with intuitive arguments, that may become 
fully rigorous) and it displays very interesting and new phenomena. 

The random energy model is defined as following. There are N Ising spins {ai, i = 1,N) that may 
take values ±1; the total number of configurations is equal to M = 2^: each configuration can be 
identified by a label s in the interval 1 — M. 

Generally speaking the Hamiltonian of the system is determined when we know all the values of 
the energies Eg, i.e. a value for each of the M configurations of the system. In many models there 
is an explicit expression for the energies as function of the configurations; on the contrary here the 
values of the Eg are random, with a probability distribution p{E) that is Gaussian: 

KE) oc exp . (8) 

The energies are random uncorrelated Gaussian variables with zero average and variance equal to N. 
The partition function is simply given by 

Zi= eM-PEs) = 2^ f pi{E)eM-m, (9) 

s=l,M ■' 

s=l,M 

The values of the partition function and of the free energy density (// = —\n{Zi)/{N0)) depend 
on the instance / of the system, i.e. by all the values of the energies Eg. It can be proved that when 
N ^ oo the dependance on / of // disappears with probability 1. As usual, the most likely value of 
// coincides with the average of //, where the average is done respect to all the possible instances of 
the system. 

In the following I will indicate with a bar the average over the instances of the system: 

Ai = j dii{I)A{I) , (10) 

where //(/) is the probability distribution of the instances of the system. Using this notation we would 
like to compute 

f = Ti (11) 

and to prove that 



(// - f? ^N^oo (12) 

The model is enough simple to be studied in great details; exact expressions can be derived also 
for finite N. In the next section we will give the main results without entering too much into the 
details of the computation. 
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2.2 Equilibrium properties of the model 

The first quantity that we can compute is Zi. We have immediately 



Zj = ^exp{-(3Es) = 2^exp(:^/?2) . (13) 



However the physically interesting quantity is ln(Z7), that is related to the average of the free energy: 



-NPf = ln{Zi) (14) 

If we make to bold assumption that 

HZj) = In (Z7) (15) 
(this approximation is usually called the annealed approximation), we find that 

I^ = iv(ln(2) + i/32) . (16) 

Generally speaking the relation eq. (15) is not justified. In this model it gives the correct results in 
the high temperature phase, but it fails in the low temperature phase. For example, as we shall see 
later, at low temperature we must have that the quantity 



is different from the physically interesting quantity 

d_ 

dp 

where {E)i is the expectation of the energy in the instance / of the problem: 



E = -4^1n {Zi) = J dn{I){E)i , (18) 



^^^^ - E.exp(-/3^,) 

In the computation of G we give to the different systems a weight equal to the partition function 
Zi that is exponentially large (or small) and strongly fluctuates from system to system. The correct 
results E is obtained using the correct (flat) weight dfi^I). 

One can argue (and it will be done later on in section 2.4) that the previous derivation gives the 
correct form for the free energy in the low (3 region. However the previous results cannot be correct 
everywhere. If we identify G with E and we compute the energy density of the model using (eq. 16), 
we find 

e = l3 (20) 



that superficially is not contradictory. However, if we use this expression for the energy to compute 
the entropy density, we find 

S{(5) = = ln(2) - ^ (21) 



This entropy density becomes negative for 

/? > /?e = ^21n(2) . (22) 
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Figure 2: The quantity ln(A''o(e))/A^ as function of the energy density e outside the interval the 
function A/o(e) is zero. 

This is a major inconsistency in a model where the configuration space is discrete and the entropy 
must be a non-negative quantity. 

The arguments we have used in deriving eq. (16) are not solid. We have computed the asymptotic 
behaviour of Z„ for positive integer n when the volume goes to infinity, but we cannot argue that the 
asymptotic behaviour is an analytic function of n. A simple counterexample is given by 

sinh(iVn) ^ exp(iV|n|) (23) 

The non-analyticity of the free energy density as functions of a parameter in the infinite volume limit 
is the essence of phase transitions. It is quite possible that at low temperatures a phase transition in 
n is present in the interval < n < 1, so that the extrapolation of the free energy density from n > 1 
to n = has nothing to do with the value of the free energy computed directly at n = 0. 

There is a simple way to find the correct result [19, 20]. The probability of finding a configuration 
of energy in the interval [E, E + dE] is given by 

M^{E)dE , (24) 

where 

M,{E) ^ 2^ exp [-^ = exp(iV(ln(2) - ^e')), (25) 

and e = E/N is the energy density. In the case of a generic system (with probability 1 when N oo) 
no configurations are present in the region where J\fo{E) << 1, i.e. for 

< = 21n(2). (26) 

Neglecting prefactors, the average total number of configurations of energy less that e (when e < 0) 
is given by 
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Figure 3: The entropy density as function of the temperature. 

pNe 1 

/ dEMo{E) ^ exp(7V(ln(2) - -e^)) . (27) 

J —oo ^ 

This number becomes exponentially small as soon e < Cc, implying that for a generic system there 
are no configurations for e < Cc- We can thus write for a generic instance of the system the following 
relation: 

p{E) « M{E) = Mo{E)e{El - E^). (28) 
The partition function can be written as 

J Af{E)exp{-PE). (29) 

Evaluating the integral with the saddle point method we find that in the high temperature region, 
i.e. for p < f3c = e~^, the internal energy density is simple given by —(3. This behaviour must end 
somewhere because we know that the energy is bounded from below by —Cc- In the low temperature 
region the integral is dominated by the boundary region E k, E^ and the energy density is exactly 
given by -Cc. 

The point where the high temperature behaviour breaks down is exactly the point where the 
entropy becomes negative. The entropy density is positive in the high temperature region, vanishes at 
f3c and remains zero in the low temperature region. In the high temperature region an exponentially 
large number of configurations contributes to the partition function, while in the low temperature 
region the probability is concentrated on a finite number of configurations [20] . 

It is interesting to notice that in the low temperature phase the free energy is higher than the 
continuation from high temperature. This is the opposite of the usual situation: in conventional mean 
field for non-disordered systems (e.g. ferromagnets) the low temperature free energy is lower than 
the continuation from high temperature ^. This correspond to a jump upward of the specific heat 
when we decrease the temperature, while in the REM (and in real glasses) there is jump downward 

^If two branches of the free energy are present usually the relevant one has the lower free energy. We shall see later 
that this is not the case in glassy systems. 
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of the specific heat. The system has less configurations of what can be inferred from the behaviour 
at entropy at high temperature: some part of the predicted phase space is missing and consequently 
the free energy is higher than the analytic continuation of the high temperature results. This entropy 
crisis is the essence of the Kauzman transition [8]. 

A more precise (and less hand-waving) computation can be done by using the representation 



HZ) = / 

Jo 



exp(-0 - exp(-^Z) ^^^^ 



and writing 

^M=tz)= E (31) 

Oil. 
,00 

The quantities Z"- can be computed exactly and in this way we can also estimate the finite correc- 
tions to the asymptotic behaviour. The details of this computation can be found in [20]. 

2.3 The properties of the low temperature phase 

It is worthwhile to study the structure of the configurations that mostly contribute to the partition 
function in the lower temperature phase. The detailed computation of the finite N corrections [20] 
shows that in the low temperature phase the average total entropy has a value that remains finite in 
the limit iV — 00: 

where the quantity m is given by 

m = pc/P = T/Tc . (33) 

The parameter m is always less the one in the whole temperature phase. 

The finiteness of the total entropy implies that the probability distribution is concentrated on 
a finite number of configurations (see eq. 45) also in the limit N ^ 00. In order to study the 
contributions of the different configurations to the free energy, it is useful to sort the configurations 
with ascending energy. We relabel the configurations and we introduce new labels such that < Ei 
for k < i. 

It is convenient to introduce the probability Wk that the system is in the configuration k; 

_ e^p{-pEk) , . 

Wk = ^ • (34) 

It is obvious that the Wk form a decreasing sequence and that 

E ^k = i- (35) 

fe=l,2^ 

Here the probability of finding an energy E^ in the interval [E,E + dE], not far from the ground state 
energy Ei is given by 

p{E) « cxp{/3^{E - E*)) = exp{f3miE - E*)). (36) 

where E* is a reference energy not far from the ground state (the parameter m controls the local 
exponential growth of the density of configurations). It is not difficult to prove [21] that, if we define 

A = /3,(Ei - E*) , (37) 
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Figure 4: The free energy density as function of the temperature and the analytic continuation of the 
high temperature result (dashed) 

the probability distribution of A is given by the Gumbel law: 

P(A) = exp(A - exp(A)). (38) 

It is somewhat more difficult 3 to prove that probability of finding one of the variables w in the interval 
[•10,10 + dw] is given 



^ ~ r(m)r(i - m) ■ ^^^^ 



where the proportionality factor can be easily found using the condition 

r-l 







dwi'{w) w = 1 . (40) 



Now it easy to compute 
As particular cases we get 



k T{k-m) 



" = r(i-m)r(fe) • 

■ , / X 2 , /■ / X s (l-m)(2-m) ,^ . 

dwu{w)w^ = 1 - m, / u{w)w^ = '-^ . (42) 

We can check that the total entropy is given by 

St = — J d'wi'{w)wln{w) , (43) 

^The proof sketched in one of the appendices. 
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/' 

Jo 



The integral 

nl 

dwu{w) (44) 

/o 

is divergent, signaling that there are an infinite number of w's. 

A detailed computation [20, 2, 22] shows that a finite number of terms dominates the sum over k. 
Indeed using the esplicite form of u{w) near w = we get that 

^ Wk = l- 0(L-^), A = ^—^ . (45) 

In this model everything is clear: in the high temperature region the number of relevant config- 
urations is infinite (as usual) and there is a transition to a low temperature region where only few 
configuration dominates. 

The fact that the low temperature probability distribution is dominated by a few configurations 
can be seen also in the following way. We introduce a distance among two configurations a and 7 as 

d{a,j) = — . (46) 

We also introduce the overlap q defined as 

?(«,7) = ^^^^^^ = l-ciV,7). (47) 
The distance squared is normalized in such a way that it spans the interval — 2. It is equal to 

• 0, if the two configuration are equal {q = 1). 

• 1, if the configuration are orthogonal {q = 0). 

• 2, iiaf = -a]{q = -l). 

We now consider the function Q{d) and P{q), i.e. the probability that two equilibrium configurations 
are at distance d or overlap q respectively. We find 



• For T > T, 



c 



Qid) = Sid-l), Piq)=S{q). (48) 

• For T < Tc 

Q{d) = {l-A)S{d-l)+AS{d), P{q) = {l-A)S{q)+A5{q-l). (49) 

where for each system (/) Aj is equal to I]fc=i,2^^fc- The average of Aj over the different 
realizations of system is equal to 1 — m. 

As soon as we enter in the low temperature region, the probability of finding two equal configura- 
tions is not zero. The transition is quite strange from the thermodynamic point of view. 

• It looks like a second order transition because there is no latent heat. It is characterized by a 
jump in the specific heat (that decreases going toward low temperatures). 
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• It looks like a first order transition. There are no divergent susceptibilities coming from above 
or from below (in short range models this result implies that there is no divergent correlation 
length). Moreover the minimum value of d among two equilibrium configurations jumps discon- 
tinuously (from 1 to 0). 

• If we consider a system composed by two replicas {a^ and a^) [23] and we write the Hamiltonian 

H{a^,a^) = H{a^) + H{a^) + Ned^{ai,a^) (50) 

for e = the free energy is equal to that of the previous model (apart a factor 2) but wc find a 
real first order thermodynamic transition, with a discontinuity in the internal energy, as soon as 
e > 0. The case e = is thus the limiting case of a bona fide first order transitions. 

These strange characteristics can be summarized by saying that the transition is of order one and a 
half, because it shares some characteristics with both the first order and the second order transitions. 

It impressive to note that the thermodynamic behavioTir of real glasses near Tc is very similar to 
the order one and a half transition of REM. We will see later that this behaviour is also typical of the 
mean field approximation to glassy systems. 

2.4 A ceireful analysis of the high temperature phase 

Let us now try to justify more carefully the correctness of eq. (15) in the high temperature phase [20]. 
We can firstly compute the quantity 

Z2 = Z] = J^^M-PEs - PE,,) (51) 

s,s' 

If we consider separately the terms with s = s' and those with s ^ s' , we easily get 

Z2 = 2^(2^ - 1) exp(2^) + 2^ exp(4^) ^ 2^^ exp(2^) + 2^ exp(4^) , (52) 
where we have used the relation 

j p{E)dE exp{mpE) = exp{m'^A) , (53) 

with 

A = (54) 

In a similar way we can compute Zf: we neglect terms that are proportional to 2~^ and we separate 
terms with all indices equal, terms with all the indices different and terms with two indices equal and 
one different; finally we find that 

= « 2=^^ exp(3A) + 3 2^^ exp((4 + 1)^) + 2^ exp(9A) . (55) 

Generally speaking, using the multifactorial formula and taking the leading contribution for each term 

Z„ = Zf = ^ ^ 2'^ eMmU) , (56) 
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where the constraint Y^\=i rrii = n\s satisfied. The previous formula can be derived by writing 

= E exp(/3 EsJ . (57) 

si---Sn a=l,n 

We now divide the n indices in / sets of size rrii and estimate the contribution from each set, that is 
given by 

Z7^=2^exp(mf^) . (58) 

The final results eq. (56) is obtained by multiplying the contribution of each set and adding the a 
prefactor (i.e. n\/{U ni=i i i^i^-)) that has a combinatorial origine. It is a simple exercise to check that 
the previous formula reduce to the one that we have written for n = 1, 2, 3. 
The previous formula can also written as 

= E E „^r' , ' exp(iV Y: G{m,)) (59) 

I {mi}''-i-i-i=l^i- i=l,l 

where 

G{mi) = ln(2) + Arnf . (60) 
For small /? and for n not too large the term I = n and rrii = l^i is dominating: it gives that 

Zn ~ exp(iVn(ln2 + A)) . (61) 

We can now use the relation 



In (Z") 



ln(Z) = hm — ^ — ^ (62) 
and apply the previous formula also for non-integer n up to n = 0. In this way one finds that 



N 2 

and therefore 



-ln2 + /l = hl2 + i/3^ (63) 



W) = -!f^-i/' (64) 



The same results can be obtained manually by writing 



This gives 



Z = 2^-PY^s + ^Ee's+--- (65) 



ln(Z) = iVln2 - - i/?^ (%^)' + + • • • (66) 



Only the first and the last term contribute in the limit N ^ oo (the energies have zero average 
and their number is 2^); a carefully evaluations shows that the remaining terms give a vanishing 
contribution. 
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2.5 The replica method 

Wc may wonder if the previous results in the low energy phase may be obtained by starting from eq. 
(56) and by partitioning the replicas using an appropriate choice of the variables that indicates 
the size of the different sets (or blocks). 

The correct choice in the high temperature region was = 1. An alternative choice (that is 
highly symmetric) is nii = m. The total number of partitions I is given by n/m. This saddle point 
gives a contribution [20, 2] 

Zn « exp(iVn(i^ + mA)) , (67) 
m 

that correspond to 

-P7 = Fim)^^ = '^ + lp^m. (68) 
m m 2 

The value of m is arbitrary. It is surprising that in the low temperature phase the correct result for 
the free energy is obtained by taking the maximum ^ of F{m) with respect to m, and the maximum 
is found in the interval — 1 (in the high temperature phase the correct result is given -F(l)). 

This prescription can looks rather strange, but we can test some consequences of this approach for 
other quantities. At this end we take two replicas of the same system and we compute the probability 
that these two configurations are the same using this approach and we compare with the known results. 

Before going on we have to make a general remark. The average over the samples of the statistical 
averages can be simply obtained by computing the average over an ensemble of n replicas of the same 
system and taking the limit n — >^ at the end. The precise statement is the following: 



{A{a)) = 



E.^(^)exp(-,jJ/(a)) 
E,exp(-/3if(a)) 



Ec.i...a„-4(c7i)exp(-/?Ea=l,ng(^a)) .gg. 
E.i....„exp(-/3Ea=l,n^^K)) 



Indeed the last term is equal to 



E. A{a) eM~f^H{a)) (E. exp(-/3iJ(a)))"-^ 
(E.exp(-/3F(a)))" 

In the limit n ^ we recover the previous equation. A similar argument tells us that 



_ln(Z)^ (71) 

where 

In (Z^) 

Here to proof can be easily done by using 1 + ?7.1n(Z) and the obvious relation 1 = 1. 

We can follow the same strategy for computing quantities that depend on a pair of replicas. Let 
us define: 

(2)/ X Esi...s„exp(-/3Ea=l,n^(Sa))'5si,s2 

'{n) = — — . (73) 

Esi...Sn ^-^Pl"/^ Ea=l,n 



*li the correct results are given by the saddle point where the replicas are partitioned in equal size blocks (with m ^ 1) 
we will say that the replica broken at one step. 
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Also in this case we find 



Indeed we can write that 



^ ^ = '"^^^^ = 1™ w^'^\n) . (74) 

fc=l,oo 



(2) . ^ _ exp{-p2E{s)) (E. expi-f32E{s))r-' 

Let us now compute W2{n) by considering the contribution of the saddle point where all the mj 
are equal to m[2]. In this case W2{n) is given by the probability of finding the replicas 1 and 2 inside 
the same block when we average over the ways in which one can divide the n replicas into n/m blocks 
of size m. The computation looks complex, but it greatly simplifies by noticing that this probability 
does not depends on the two replicas (1 and 2) that we are actually considering: it coincides with 
the probability that for a given partition, two random replicas stay in the same block. A simple 
computation tell us that this probability is given by 

n(m - 1) 

Indeed the first replica is arbitrary and the second replica is any of the remaining m — 1 replicas of 
the block. The denominator is the total number of unordered pairs of replicas. In the limit n — >^ we 
recover the correct result: 

W2 = 1 — m . (77) 

In the similar way we find 

, n(m - l)(i-n - 2) l- ;n)(2-m) r(3 - /») 

W3 = lim = '-^ ^ = —— (78) 

n-*o n{n-l){n-2) 2 r(l - m)r(3) ^ ^ 

The computation can be trivially done for any k and from the moments we recover the distribution 
probability ^{w). 

At this stage it is not clear why we have chosen this form of the nii and not a different one. 
However it is clear from the previous computations that the choice of the rrii codes the probability 
distribution of the w and the situation where the logarithms of the w^s are uncorrelated identically 
distributed variables with an exponential distribution correspond to nii = m. Different situations may 
correspond to different forms of the nii, but the one that is relevant here corresponds to constants 
rrii. For other purposes, i.e. for finding finite size corrections or exponentially vanishing probabilities, 
different choices of the nii are relevant [24] . It may be interesting to remark that the one that we have 
taken here is in some sense the maximally symmetric one if we esclude the trivial choice mj = 1. 



2.6 Dynamical properties of the model 

The dynamical properties of the model can be easily investigated in a qualitative way. Interesting 
behavior is present only in the region where the value of N is large with respect to the time, but it 
will not be discussed here. 

We consider here a single spin flip dynamics. More precisely we assume that in a microscopic time 
scale scale the system explore all the configurations that differ from the original one by a single spin 
flip and it moves to one of them (or remain in the original one) with probability that is proportional 
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to exp(— This behaviour is typical of many dynamical process, like Glauber dynamics, Monte 
Carlo, heath bath. 

In this dynamical process each configuration C has N nearby configurations C to explore. The 
energies of the configurations C are uncorrelated to the energy of C, so that they are of order A^^/^ in 
most of the case. A simple computation gives that lowest energy of the configurations C would be of 
order — (A'"ln(Ar))^/^. The important point is that the corresponding energy density (-(ln(Ar)/Ar)V2) 
vanishes in the large N limit. 

If the configuration C has an energy density e less that zero, the time needed to do a transition 
to a nearby configuration will be, with probability one, exponentially large. For finite times at large 
N a configuration of energy e < is completely frozen. Only at times larger that exp(/3eA'^) it may 
jump to a typical configuration of energy zero. At later times different behaviours arc possible: the 
configuration comes back to the original configuration of of energy e or, after some wandering in the 
region of configurations of energy density ~ 0, it fells in an other deep minimum of energy e' . A 
computation of the probabilities for these different possibilities has not yet been done, although it 
should not too difficult. 

The conclusions of this analysis are quite simple. 

• Every configuration of energy e < is a deep local minimum of the Hamiltonian: if we flip one 
spin, the energy increases of a large quantity eN). 

• If we start from a random configuration, after a time that is finite when N — ^ oo, the system 
goes to a configuration whose energy is of order — ln(A'")^/^ and stops there. 

• If we start from a random configuration, only at exponentially large times the system will reach 
an energy density that is different from zero. 

3 Models with partially correlated energy 
3.1 The definition of the models 

The random energy model (REM) is rather unrealistic in that it predicts that the energy is completely 
upset by a single spin flip. This feature can be eliminated by considering more refined models, e.g. 
the so called p-spins models [25, 26], where the energies of nearby configurations are also nearby. We 
could say that energy density (as function of the configurations) is not a continuous function in the 
REM, while it is continuous in the p-spins models, in the topology induced by the distance, (46). In 
this new case some of the essential properties of the REM are valid, but new features are present. 

The Hamiltonian of the p-spins models depends on some control variables J, that have a Gaussian 
distribution and play the same role of the random energies of the REM and by the spin variables a. 
For p = 1, 2, 3 the Hamiltonian is respectively given by 

H'Act) = Yl Ji^i ' (79) 
i=l,N 

I 

Hj{T) = Ji,k(^i(^k , (80) 

i,k=l,N 
I 

i,k,l=l,N 
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where the primed sum indicates that all the indices are different and ordered. The variables J must 
have a variance of (7V(i-p)/2) if the system has a non-trivial thermodynamical limit. Here we will 
study the hard spin models [7] where crj = ±1 

It is possible to prove by an explicit computation that, if we send first N ^ oo and later p oo, 
one recover the REM [25]. Indeed the energy is normalized in such a way that it remains finite when 
p — ^ oo; however the differences in energy corresponding to one spin flip are of order ^/p for large 
p (they ar order N in the REM), so that in the limit p — ^ oo the energies in nearby configurations 
become uncorrelated and the REM is recovered. 



3.2 The replica solution 

The main new property of the p spin model is the presence of a correlation among the energies of nearby 
configurations. This fact implies that if C is a typical equilibrium configuration, all the configurations 
that differ from it by a finite number of spin flips will differs energy by a bounded amount, also for 
very large N. The equilibrium configurations are no more isolated (as in REM), but they belongs to 
valleys: the entropy restricted to a single valley is proportional to N and it is an extensive quantity. 

We can now proceed as before ^. The computation of Z„ = can be done by introducing n 
replicas [2, 26]. After the Gaussian integration over the J one gets 

Z„ = ^exp[liV/32 ^ Qa,b{^r] , (81) 

S y a,b=l,n I 

where S denotes the set of all nN a variables and 

Qa,m = AT-i J2 ('a{i)(Tb(i) . (82) 
i=l,N 

If we introduce the Lagrange multiplier ^ we find that the previous expression reduces to 



= JdqdAj^ewllNp^ ^ g^, + iV Yl KAQa,b{^)-qa,b)] , (83) 

E y a,b=l,n a,b=l,n ) 



where both q and A are ra x n symmetric matrices and the integrals dq and dA runs over these matrices 
The sum over the spins can be done in each point independently from the other point using the 
relation 

exp(Ar A„,ft(g„,ft(E)) = exp j ^ ^ K,^a a{i)cj^{k) \ . (84) 



Finally we get the result: 



/ 



dgexp I I Y €,b + - Qa,b) I I (85) 

a,b=l,n 



different a well studied model is the spherical model [11, 9, 23], that has the same Hamiltonian: however the spins 
arc real variables that satisfy the constraints f ; = A*'. The spherical model has many features in common with the 
hard spin model and sonic computations arc simpler. 

®In this case we must use some technical tools that are also used in the study of the infinite range ferromagnetic 
model, described in the appendix. 

^The integrals over the variables Aa^b is along the imaginary axis and goes from —ioo to +ioo. 
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where 

exp(G(A)) = ^exp( ^ K,b(yaOb) (86) 

f a,b=l,n 

In the hmit ^ oo the previous integral is dominated by the saddle point. We have to find out the 
solution of the equations 



^a,b , 



2 

qa,b = (o-aO-fo)A , (87) 

where 

, , dG{A) E<7CXp(Ea,6=l,n^a,bf^aCTb)0-cCrd 

dAc,d Eaexp(Ea,6=l,n^a,60-aO-6) 

The free energy density (per replica) is given by 

fn = ih\\p" E <6+G(A)- E K^iaA ■ m 

\ a,b=l,n a,b=l,n J 



The internal energy density is given by deriving Pfn{P) with respect to /?. One finds: 

^ = ^/3 E €,b (90) 

a,b=l,n 

In the high temperature phase we have that the only solution to the saddle equations is 

qa,b = for a / 6 
Qafi = 1 for a = 6 . (91) 

The corresponding value of the internal energy is given by 

E = -P . (92) 

Also in this model the high temperature results cannot be true for all (3 because it leads to a negative 
entropy at low temperatures. 

dm dq ^ ^^^"^ 

Following the approach of the previous section we divide the n indices into sets of m indices and 
wc put Qa^b = 9 if a and b belongs to the same set and = if a and b do not belong to the same 
set. (Obviously qa,a must be equal to sz^ = 1)- In the same way A^^b = A, if a and b belongs to the 
same set, and Aa,b = 0, if a and b do not belong to the same set. Now the computation in eq. (88) 
can be done independently for each set of m replicas. The typical quantity we have to evaluate is 

E^^P E \^(^a(^b = E^^P E 

\a,b=l,m ) " \ \a=\,m 

= I dn{h) exp I VXhi E I = / Mh) cosh( VA/i)"* (94) 

\ a=l,m J 
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Figure 5: The qualitative behaviour quantity qea and m (dashed hne) as functions of the temperature 
(in units of Tc) in the p spin model for p > 2. 



where h is a Gaussian variable with variance one: 



dti{h) = (27r)-i/2 exp(-/iV2) 



(95) 



We finally find the equations 



A 



/ diiih) cosh(\/A/i)"^) tanh2(^/A/i) 
/ d/i(/i) cosh(\/A/i)'" 



(96) 



For each value of m we can find a solutions (or more solutions) to the previous equations and we 
can compute the corresponding free energy. If we use the same prescription than in the REM the 
correct solution is the maximum of the free energy as function of m for < m < 1. In other words 
we obtain a function F{m) that we have to maximize. It may looks strange to maximize the free 
energy, but this is what we need if the final free energy must be higher of the one computed in the 
high temperature phase, i.e. F{1). Moreover there is a rigourous theorem, recently proved by Guerra 
[27] , that the true free energy / satisfies the relation 



f > max F(m) 

~ 0<m<l ^ ' 



(97) 



so maximization^ and not minimization^ of the free energy is the natural prescription in this situation 
[2]. 

In general, after having eliminated A the free energy can be written as function of m and q {f{q, m)) 
and the value of these two parameters satisfy the stationarity equations: 
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Figure 6: The quantity qea and m (dashed line) as functions of the temperature in the p spin model 
(in units of Tc) for p = 2 (i.e. the SK model). 

If we consider the case p > 2, we finds that for T < there is a solution where the value of q 
at the maximum of f{q,m) (that we will call qea) is different from zero. The value of qea jumps at 
Tc to a non-zero value. For large values of p the quantity qea {qEA would be 1 in the REM) is of 
order 1 — exp(—APp), while the parameter m has the same dependence on the temperature as in the 
REM, i.e. it is equal to 1 at the critical temperature and has a linear behaviour al low temperature. 
When p is finite m is no more strictly linear as function of the temperature. The thermodynamical 
properties of the model are the same as is the REM: a discontinuity in the specific heat, with no 
divergent susceptibilities. 

The case p = 2 (the Sherrington Kirkpatrick model [28] that is relevant for spin glasses) has a 
different behaviour especially near the critical temperature (see fig. 6). 

3.3 The physical interpretation 

Also in these models for each instance / we can define the function Pi{q) and compute its average 
over the different realization of the system. The computation goes on exactly in the same way as in 
the REM. For example let us consider 



where the values of the replica indices (i.e. a and b) do not matter, as far as they are different. In the 
same way as before we average over all the ways in which we can assign the two replicas to different 
blocks. We find 




(98) 




(99) 



where qa,b 



is the average of qa^b over the n! permutations of the replicas. 
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In the same way 

= ) = <^ = (1 - ^khA ■ (100) 

Putting everything together one finds for the average over the different instances 

PAm) = P{q) = m5{q) + (1 - m)5{q - qea) ■ (101) 
One can painfully reconstruct all the probability distributions of the w' shy studying quantities like 



Pi{q)<l) =<llAY.'^k=W^ (102) 
^ k 

and doing the appropriate combinatorial estimates, at the ends one finds (exactly in the same way as 
in the REM) that they are proportional to exp(— /3F/j), where the are distributed according to the 
distribution eq. (36) . 

In order to develop a formalism useful to discuss the physical meaning of these results it is conve- 
nient to introduce the concept of pure states in a finite volume^!]. This concept is crystal clear from 
a physical point of view, however it can be difficult to state it in a rigorous way (i.e to prove existence 
theorems). 

We consider a system in a box of linear size L, containing a total of N spins. We partition 
the configuration space in regions, labeled by a, and we define averages restricted to these regions 
[29, 30, 4]: these regions will correspond to our finite volume pure states or phases. It is clear 
that in order to produce something useful we have to impose sensible constraints on the form of these 
partitions. We require that the restricted averages on these regions are such that connected correlation 
functions are small at large distance x in a short range model (or when the points are different in an 
infinite range model). This condition is equivalent to the statement that the fiuctuation of intensive 
quantities ^ vanishes in the infinite volume limit inside a given phase. 

In a ferromagnet the two regions are defined by considering the sign of the total magnetization. 
One region includes configurations with a positive total magnetization, the second selects negative 
total magnetization. There are ambiguities for those configurations that have exactly zero total 
magnetization, but the probability that such a configuration can occur is exponentially small at low 
temperature. 

In order to present an interpretation of the results we assume that such decomposition exists 
also each instance of our problem. Therefore the finite volume Boltzmann-Gibbs measure can be 
decomposed in a sum of such finite volume pure states. The states of the system are labeled by a: we 
can write 




a 

with the normalization condition 

Y.Wa = l. (104) 
a 

The function Pj{q) for a particular sample is given by 

^Intensive quantities are defined in general as 6 = 'Yla-i ,where the functions Bi depend only on the value of CTj 
or from the value of the nearby spins. 
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Pji^l) = WaWfj5{<ia,[3 - q) , (105) 

where Qq,^^ is the overlap among two generic configurations in the states a and p. 

Given two spin configurations {a and r) we can introduce a natural concept of distance by 

1 ^ 

d\a,T)^-Y.(<',-nf , (106) 
1=1 

that belongs to the interval [0 - 1], and is zero only if the two configurations are equal. In the 

thermodynamical limit, i.e. for N oo, the distance of two configurations is zero if the number of 
different spins remains finite. The percentage of different a's, not the absolute number, is relevant in 
this definition of the distance. It is also important to notice that at a given temperature P~^, when 
N goes to infinity the number of configurations inside a state is extremely large: it is proportional to 
exp{NS{f3)), where S{P) is the entropy density of the system). 

We expect that finite volume pure states will enjoy the following properties that likely characterizes 
the finite volume pure states: 

• When N is large each state includes an exponentially large number of configurations^ 

• The distance of two different generic configurations Ca and (the first belonging to state a and 
the second to state 7) does not depend on the Ca and C^, but only on a and /3. The distance 
da,i3 among the states a and /3, is the distance among two generic configurations in these two 
states. The reader should notice that with this definition the distance of a state with itself is 
not zero. If we want we can define an alternative distance: 

Da,l3 = da,l3 - ^ {da,a + d/?,/?) , (107) 

in such a way that the distance of a state with itself is zero (-Da,a = 0). 

• The distance between two configurations belonging to the same state a is strictly smaller than 
the distance between one configuration belonging to state a and a second configuration belonging 
to a different state /?. This last property can be written as 



da,a < da,l3 ■ (108) 

This property forbids to have different states such that -Da,/3 = 0, and it is crucial in avoiding 
the possibility of doing a too fine classification 

• The classification into states is the finest one that satisfies the three former properties. 



The first three conditions forbid a too fine classification, while the last condition forbids a too 
coarse classification. 

®We warn the reader that in the case of a glassy system it is not possible to consider N ^ oo limit of a given finite 
volume pure state: there could be no one to one correspondence of states at N and those at 2A''.. 

^°For example if in a ferromagnet at high temperature we would classify the configurations into two states that we 
denote by e and o, depending on if the total number of positive spins is even or odd, we would have that de,e = de,o = do,o- 



25 



For a given class of systems the classification into states depends on the temperature of the system. 
In some case it can be rigorously proven that the classification into states is possible and unique 
[31, 32, 33] (in these cases all the procedures we will discuss lead to the same result). In usual 
situations in Statistical Mechanics the classification in phases is not very rich. For usual materials, in 
the generic case, there is only one phase. In slightly more interesting cases there may be two states. 
For example, if we consider the configurations of a large number of water molecules at zero degrees, we 
can classify them as water or ice: here there are two states. In slightly more complex cases, if we tune 
carefully a few external parameters like the pressure or the magnetic field, we may have coexistence 
of three or more phases (a tricritical or multicritical point). 

In all these cases the classification is simple and the number of states is small. On the contrary in 
the mean field approach to glassy systems the number of states is very large (it goes to infinity with 
N), and a very interesting nested classification of states is possible. We note "en passant" that this 
behavior implies that the Gibbs rule is not valid for spin glasses. 

The results we have obtained on the probability distribution of the overlap can be interpreted 
by saying that the system has many equilibrium states with weight that have the same probability 
distribution as a REM with the appropriate value of m such that the average value of the local 
magnetization squared is inside each state is just qea- 

E»=i,jv(?-^a(0)^ nno^ 
-j^ = QEA (109) 

Different states have magnetization that points toward orthogonal directions: 



N 



(110) 



if a 7^ 7. 

The only difference with the REM is that qea is less than 1. The states in the REM are composed 
by a single configuration, in the p spin models local fluctuation are allowed and the states contain an 
exponentially number of conflgurations (local fluctuation are allowed) and have a non-zero entropy. 



3.4 The two susceptibilities 

Let us consider a model where the Hamiltonian is of the form 

Ho{a)-hY^ai (111) 

i 

If the variables are interpreted as spins, h is the magnetic field. A crucial question is what happens 
when we change the magnetic field. 

In order to simplify the analysis let us suppose that we stay in a system where, after we average 
on the disorder, ioi i ^ k: 

{ciOk) = = . (112) 

This property is valid in the p-spin models and also in simple models of spin glasses: e. g. in the 
Edwards Anderson model [1] where 

Ho{(t) = ^ Ji,k<7i(Tk , (113) 

i,k 

^^The Gibbs rule states that in order to have coexistence of n phases (n-critical point), we must tune n parameters. 
Here no parameters are tuned and the number of coexisting phases is infinite! 
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where the sum over i and k runs over the nearest neighbours and the J are random variables with 
even probability distribution. 

A simple computation [2] shows in general that the average equilibrium magnetic susceptibility is 
just given by 

Xeg = P J dq P{q)il - q) . (114) 

Indeed 

i,k=l,N 

The terms with i ^ k do not contribute after the average over the systems (as consequence of eq. 
112): the only contribution comes from the terms where i = k. We finally obtain 

NXeq = E (1 - ii^i) f = (l-E^«a,7) • (116) 

i=l,N \ a,y / 

In p spin models using the form of the function P{q) or using directly the properties of the w^s and of 
the g's we find 

Xeq = /3(1 - (1 - rnjqEA) = /3(1 - qEA + mqEA) (117) 
It is interesting to note that we can also write 

Xeg = P{1 - QEA) + ^Y^^aW^iMa - M^f (118) 

where Mq, is the total magnetization is the state a 

Ma = Y.^a(i))a (119) 

i 

The first term (i.e. /3(1 — qEA)) has a very simple interpretation: it is the susceptibility if we 
restrict the average inside one state and it can be this identified with xlr- The send term gives Xirr'- 

Xirr = PY '^a'Wj{qEA - qa,j) = PmqEA , (120) 
a,7 

where the last equality is valid in the p-spin models. In these models (for p > 2) Xirr jumps at 
the critical temperature (x remains continuos): in spin glasses (p = 2) qEA vanishes a the critical 
temperature and both x and Xirr are continuos in agreement with the experimental results (fig. 7). 

The physical origine of Xirr is clear. When we increase the magnetic field, the states with higher 
magnetization become more likely than the states with lower magnetization: this effect contributes 
to the increase in the magnetization. However the time to jump to a state to an other state is very 
high (it is strictly infinite in the infinite volume limit and for infinitesimal magnetic fields where we 
can neglect non-linear effects are neglected): consequently the time scales relevant for xlr and Xeq 
are widely separated. 

If we look to real systems (e.g. spin glasses) both susceptibilities are experimentally observable. 

• The first susceptibly {xlr) is what we measure if we add an very small magnetic field at low 
temperatures. The field should be small enough in order to neglect non-linear effects. In this 
situation, when we change the magnetic field, the system remains inside a given state and it is 
not forced to jump from a state to an other state and we measure the ZFC (zero field cooled) 
susceptibility, that corresponds to xlr- 
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Figure 7: The experimental results for the FC (field cooled) and the ZFC (zero field cooled) magneti- 
sation (higher and lower curve respectively) vs. temperature in a spin glass sample (Cu87Mni3.5) for 
a very small value of the magnetic field H =1 Oe (taken from [34]). For a such a low field non-linear 
effects can be neglected and the magnetization is proportional to the susceptibility. 

• The second susceptibility (Xeq) can be approximately measured by doing cooling the system 
in presence of a small field: in this case the system has the ability to chose the state that is 
most appropriate in presence of the applied field. This susceptibility, the so called FC (field 
cooled) susceptibility is nearly independent from the temperature (and on the cooling rate ^^) 
and corresponds to Xeq- 

Therefore one can identify XLR and Xeq with the ZFC susceptibility and with the FC susceptibility 
respectively. The experimental plot of the two susceptibilities is shown in fig. (7). They are clearly 
equal in the high temperature phase while they differ in the low temperature phase. 

The difference among the two susceptibilities is a crucial signature of replica symmetry breaking 
and, as far as I known, can explained only in this framework. A small change in the magnetic field 
pushes the system in a slightly metastable state, that may decay only with a very long time scale. 
This may happens only if there are many states that differs one from the other by a very small amount 
in free energy. 

3.5 The cavity method 

The cavity method is a direct approach that in principle can be used to derive in an esplicite way all 
the results that have been obtained with the replica method^^. The basic idea is simple: wc consider 
a system with N spins {i = l,N) and we construct a new system with N + 1 spins by adding an extra 
spin (at i = 0). We impose the consistency condition that the average properties of the new spin are 
the same of that of the old ones [2, 35]. 

^■^The nearly independence of the field cooled magnetization on the cooling rate can be used to argue that field cooled 
magnetization is near to the equilibrium one; on the contrary, if the field cooled magnetization would have been strongly 
dependent on the cooling rate, the statement that it correspond to the equilibrium magnetization would be quite doubtful 

^^Not all the results have been actually derived. 
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In order to lighten the notation I will write down the formulae only for p = 2, but it easy any to 
extend the computation to other models. The Hamiltonian of the new spin is; 



i=l,N 

If we suppose that the spins (jj have vanishing correlations and we take care that each individual term 
is small, we find that 

uiq = (do) = tanh(/3/i) . 

h= Yl '^o,imi , (122) 
i=i,Ar 

where rrii denotes the magnetization of the spin ai before we add the spin 0. When the variables J 
are random (or the variables rrii are random), the central limit theorem implies that /t is a Gaussian 
random variable with variance ^ 

= QEA = • (123) 

If we impose the condition that the average magnetization squared of the new point is equal to that 
of the old points, we arrive to the consistency equation: 

qEA = ^Q = j d^q^^{h)i3W{l3h) (124) 

where diJ,q^^{h) denotes a normalized Gaussian distribution with variance qea- It is easy to check 
that the increase in the total free energy of the system is 

AF{h) ^ -HcosHPh)) ^^25) 

In this way we have derived the replica symmetric solution that corresponds to m = 0. 

If replica symmetry is broken the spins are uncorrelated within one states, but if we do not separate 
the states, we find strong correlations. 

The correct computation goes as follows. We suppose that in the system with N spins we have 
a population of states whose total free energies are distributed (when is not far from a given 
reference value F* that for lighten the notation we take equal to zero) as 

A/'jv(i^jv) oc exp(^mFjv) (126) 

When we add the new spin, we will find a value of the field h that depends on the state a. We can 
now consider the conjoint probability distribution of the new free energy and of the magnetic field. 
We obtain 

J\fN+i{F, h) = j dhPs{h) J dFNJ\fNiFN)S{F -Fn- AF{h)) (127) 

where Ps{h) is the probability distribution of the effective magnetic field produced on the spin at 0, 
for a generic state and it is still given by duq^^^h). It is crucial to take into account that the new 
free energy will differs from the old free energy by an energy shift that is h dependent. If we integrate 
over Fjv and we use the esplicite exponential form for A/jv(FAr) we find that 

MN+iiF, h) oc exp(/3mF) j dhPs{h) exp{pmAF{h))) oc exp(/3mF)Piv+i(/i) (128) 
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The probability distribution of the field at fixed value of the free energy is given by 

Pjv+i(/i) cx Ps{h) exp(/3mAF(/i))) = dfiq^^{h) cosh(/?/i)'" . (129) 

and it obviously different from Ps{h) as soon as m 7^ 0. In this way we find the consistency equation 
of the replica approach for qea- 
A few comments are in order: 

• The probability distribution of h at fixed value of the free energy of the N spins system {Ps{h)) 
is not the probability distribution of h at fixed value of the free energy of the N + 1 spins system 
PN+i{h): the two free energies differs by an h dependent addictive factor and they do not have a 
flat distribution (as soon as m 7^ 0). The probability distribution of h at a fixed value of the free 
energy of the N spins system is Gaussian, but the probability distribution of h at fixed value of 
the free energy of the + 1 spins system is not a Gaussian. 

• Only in the case were A/7v(-^iv) is an exponential distribution MN+i{P,h) factorizes into the 
product of an F and an h dependent factor and the A/iv+i(-F) has the same form of Mn{F). 
Self-consistency can be reached only in the case of an exponential distribution for A/Ar(-^iv)- 

• The equations do not fix the value of m. This in natural because (as we shall see later) we can 
write them also in the case where the free energies densities we consider are different from that 

of the ground state {F* is not near to the ground state). In this case they do correspond to 
the distribution of the free energies inside metastable states that arc characterized by a different 
value of m than the ground state. These equation will be useful for the computation of the 
complexity. 

It is appropriate to add a last comment. The computation we have presented relates the magne- 
tization a spin of the systems with N + 1 spins to the magnetizations of the system with N spins: 
they are not a closed set of equations for a given system. However we can also write the expression of 
magnetization at zero as function of the magnetizations of the system with + 1 spins, by computing 
the variations in the magnetization in a perturbative way. Let us consider the case p = 2 and let us 
denote by m the magnetization of the old system (iV spins) and by m' the magnetization of the new 
system (AT -|- 1 spins). Perturbation theory tell us that 

m[^mi + Jo^irriQ^^ = rui + Jo,i"io/3(l - (m-)^ . (130) 

Using the previous formula we get the TAP equations [36, 2]: 

rriQ = tanh(/3/i) 

h= '^o,imi « '^o,imi -m'oY, •^o,i/3(l - "^•) (131) 

i=l,N i=l,N i 

(132) 

sumi=i^NJo,im'i - mo/3(l - qea) 

where {N + l)qEA) = J2i=o,N we have used the fact that Jq j = N~^. A detailed computation 

show that the free energy corresponding to a solution of the TAP equations is given by the TAP free 
energy. 

F[m] = Y Ji,kmimk - iV/3(l -qf-T^ S{mi) , (133) 

i<k i 



30 



where S{m) is the usual single spin entropy: 

1 + m 1 + m 1-m ,1 - m 
-Sim) = ) + ) • (134) 

It is important to note that the solutions of the TAP equations are also stationary 
free energy: using the relation 

dqEA _ m 
drrii N 

the TAP equations can be written as 

dF[m] ^ 
drrii 

4 Complexity 

It would be interesting to characterize better the free energy landscape of the models described in the 
previous section, especially in order to understand the dynamics. Indeed we have already seen that 
in the REM the system could be trapped in metastable configurations. In models where the energies 
are correlated the situation is more complicated; moreover in realistic finite dimensional models there 
are still further subtleties. 

Although the word metastable configuration has a strong intuitive appeal, we must define what 
a metastable configuration is in a more precise way. There are two different (hopefully equivalent) 
definitions of a metastable state or valley: 

• From an equilibrium point of view a valley is a region of configuration space separated by the 
rest of the configuration space by free energy barriers that diverge when iV — > oo. More precisely 
the system, in order to go outside a valley by moving one spin (or one particle) at once, must 
cross a region where the free energy is higher than that of the valley by a factor that goes to 
infinity with A^. 

• From the dynamic point of view a valley is a region of configuration space where the system 
remains for a time that goes to infinity with N. 

The rationale for assuming that the two definitions are equivalent is the following. We expect that for 
any reasonable dynamics where the system evolves in a continuous way (i.e. one spin fiip at time), the 
system must cross a configuration of higher free energy when it goes from a valley to an other valley. 
The time for escaping from a valley is given by 

T ~ To exp(/?AF) (137) 

where AF is the free energy barrier 

It is crucial to realize that in infinite range models valleys may have a free energy density higher 
that that of equilibrium states. This phenomenon is definitely not present in short range models. Two 
equilibrium states with infinite mean life must have the same free energy. 

The proof of this statement is simple. Indeed let us suppose that the system may stay in two 
phases (or valleys) that we denote as A and B. If the free energy density of B is higher than that of A, 
the system can go from B to A in a continuos way, by forming a bubble of radius R of phase A inside 

^''in kinetically constrained models, where some local movements are forbidden, we can have dynamical valleys that 
do not correspond to valley from the equilibrium point of view. 



points of the TAP 
(135) 

(136) 
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phase B . For example, if we take a mixture of H2 and O2 at room temperature, the probability of 
a spontaneous temperature fluctuation in a small region, that leads to later ignition and eventually to 
the explosion of the whole sample, is greater than zero (albeit quite a small number), and obviously 
it does not go to zero when the volume goes to infinity. This argument does not work in mean field 
models where in some sense surface effects are as important as volume effects (when D is large R^, 
the volume, is not so different from R^~^, the surface). 

We have two possibilities open in positioning the predictions of mean field theory concerning the 
existence of real metastable states: 

• We consider the presence of these metastable state with infinite mean life an artefact of the 
mean field approximation and we do not pay any attention to them. 

• We notice that in the real systems there are metastable states with very large (e.g. much greater 
than one year) mean life. We consider the infinite time metastable states of the mean field 
approximation as precursors of these finite mean life states. We hope (with reasons) that the 
corrections to the mean field approximation will give a finite (but large) mean life to these states 
(how this can happen will be discussed in the next section). 

In these notes I will explore the second possibility, that seems to be much more fruitful than the first 
one. 

4.1 The basic definitions 

Before discussing the difficulties related to the definition of the complexity in short range models, we 
must see the main definitions that are correct in the mean field approach. 

The basic ideas are quite simple [11, 9, 10, 23, 37, 38, 39, 40]. In principle we proceed in a way 
similar to the construction of equilibrium states described in section 3.3: we partition the whole 
configuration space into valleys. If we call the contribution of each valley to the partition function, 
the corresponding free energy is given by 

Z„ = exp(-/3F«) (138) 

This definition does not give us a practical way to find the valleys. An alternative approach, that 
should be hopefully equivalent, is the following. In many case one can prove that the magnetization in 
a given valley should satisfy some equations, e.g. the TAP equations of the previous sections: valleys 
may be identified with solutions of the TAP equations and their free energy is given by the TAP free 
energy^^.. Generally speaking in a system of N spins we can introduce a free energy functional F[m\ 
that depends on the local magnetizations m{i) and on the temperature. Only in the mean field case 
F[m] is given by the TAP free energy. 

We suppose that at sufficiently low temperature the functional F[m] has many local minima (i.e. 
the number of minima goes to infinity with the number (A'^) of spins). Exactly at zero temperature 
these local minima coincide with the local mimima of the potential energy as function of the coordinates 
of the particles. Let us label then by an index a. To each of them we can associate a free energy Fa 

^^If the surface tension among phase A and B is finite, has happens in any short range model, for large R the volume 
term will dominate the free energy difference among the pure phase B and phase B with a bubble of A of radius R. This 
difference is thus negative at large R, it maximum will thus be finite. A finite amount of free energy in needed in order 
to form a seed of phase A where the spontaneous formation of phase A will start. 

^^At zero temperature one could try to identify valleys with the minima of the Hamiltonian, that are called inherent 
structures in the glass community [6]. 
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and a free energy density = F^/N . In this way the vaUeys are associated to local minima of the 
free energy functional. 

In this low temperature region we suppose that the total free energy of the system can be well 
approximated by the sum of the contributions to the free energy of each particular local minimum. 
We thus find: 

Z = eM-PNfs) = J2 exp(- Wa) • (139) 

a 

When the number of minima is very high, it is convenient to introduce the function M{f,T,N), 
i.e. the density of minima whose free energy is near to /. With this notation we can write the previous 
formula as 

Z = l dfeM-PNf)M{f,T,N). (140) 

In the region where TV is exponentially large we can write 

Mif, T, N) ^ exp(iVE(/, T)), (141) 

where the function E is called the complexity or the configurational entropy (it is the contribution to 
the entropy coming from the existence of an exponentially large number of locally stable configura- 
tions) . 

The minimum (maximum) possible value of the free energy is given by fm{T) (/M(r)).The relation 
(141) is valid in the region fm{T) < f < fniT). Outside this region we have that M{f,T) = 0. It ah 
cases known T,{fm(T),T) = 0, and the function E is continuous at fm- On the contrary in mean field 
models it happens frequently that the function E is discontinuous at fm 

For large values of we can write 

//m 
d/exp(-Ar(/3/-E(/,r)). (142) 

We can thus use the saddle point method and approximate the integral with the integrand evaluated 
at its maximum. We find that 

Pfs = mm $(/) = Pf* - T), (143) 

where the potential $(/) (that will play a crucial role in this approach) is given by 

$(/) = /3/-E(/,T). (144) 

(This formula is quite similar to the well known homologous formula for the free energy, i.e. /?/ = 
miuEiPE — S{E)), where S{E) is the entropy density as function of the energy density.) 
If we call /* the value of / that minimize $(/). we have two possibilities: 

• The minimum /* is inside the interval and it can be found as solution of the equation /? = dH/df. 
In this case we have 

= pf* - E*, E* = E(/*, T). (145) 

The system may stay in one of the exponentially large number of possible minima. The number 
of minima where is convenient for the system to stay is exp(A^E*) . The entropy of the system is 
thus the sum of the entropy of a typical minimum and of E*, i.e. the contribution to the entropy 
coming from the exponential large number of microscopical configurations. 
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• The minimum is at the extreme value of the range of variabihty of /. We have that /* = fm and 
$ = /m- In this case the contribution of the complexity to the free energy is zero. The different 

states that contribute the free energy have a difference in free energy density that is of order 
N^'^ (a difference in total free energy of order 1). Sometimes we indicate the fact that the free 
energy is dominated by a few different minima by say the replica symmetry is spontaneously 
broken [2, 22]. 

From this point of view the behaviour of the system will crucially depend on the free energy 
landscape [41], i.d. the function S(/, T), the distance among the minima, the height of the barriers 
among them... 

4.2 Computing the complexity 

We have seen that complexity counts the number of metastable states. It would be interesting to 
compute the complexity in a direct way without having to count all the metastable states (an impos- 
sible task for large N) An interesting route to the evaluation of complexity consists in introducing 
new artificial couplings and to consider the behaviour of the systems in these conditions. This new 
approach works also in cases where the free energy functional is not known in an exact way, so that 
its minima cannot be computed. 

The basic idea is to start from an equilibrium configuration and to explore the configuration space 
phase around it [38, 39]. If we can define in some way the entropy (Sy) of the valley around a given 
equilibrium configuration, we have that 

S = ^* + S^ , (146) 

where S is the total entropy of the system and S* is the equilibrium complexity. 

More precisely, we study a system of N interacting variables ai, i = 1,...,N, with Hamiltonian 
H{a), q{a, r) is an overlap function, i.e. an intensive measure of similarity among the configurations a 
and r. We can consider a reference equilibrium configuration a, that produce a fixed external potential 
on a replica r. The partition function of the second system is: 

Z{a, e) = ^ exp [-PH{t) + PeNq{T, a)] , (147) 

T 

and 

r{e) = -{N/3)-\logZ{a,e))^ (148) 

is the e-dependent free energy {{■)^ denotes the average over the variables a). The new term in the 
Hamiltonian, for eps sufficiently large, forces the variables r to be near the variables a and produces 
a quenched disorder for the variables r. By changing the value of e we can to explore the phase space 
around a given equilibrium configuration a. At the end we average over a the logarithm of the a 
dependent free energy. 

The quantity r(e) is well defined and it may be computed in also in numerical simulations. However 
it is interesting to evaluate it in mean field models, where analytic computations are possible. The 
analytic computation of r(e) can be done by considering l+s replicas: the corresponding Hamiltonian 
is 

Hs{a) = H{ai) + ^ H{aa) + eiV E '?(^'^«) (1^9) 

a=2,s a=2,l+s 
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Figure 8: Different sliapes of the function W for various temperatures: the higher curves correspond 
to higher temperatures. 



Here ai plays the role of a and the da (for a = 2, 1 + s) are s replicas of the r variables. The quenched 
hmit (where there is no feedback reaction of the r variables on the a variables) is obtained in the hmit 

r(e) = lim ^In ( [] Y.^M-f3Hs{a))] (150) 



In models where we have to perform the average over the instances of the problems, we have to rephcate 
n times the s-rephcated system: we have to take n x (1 + s) variables and to study simultaneously the 
limit s — > and n ^ 0. Fortunately enough, in the high temperature phase, where we are interested 
to the computation of the complexity, we do not need to break the replica symmetry (at least for the 
computation of the complexity) and we obtain the correct results already for n = 1. 

In the following we will study the phase diagram of the model in the e — T plane. Explicit 
computations can be found in the literature mainly for the p-spins spherical model, but the conclusions 
have a general validity [38, 39]. At this end it is convenient to define the Legendre transform of r(e), 
defined as 



W{q)=T{e{q))-e{q) 
dW{q) 



dq 



eiq) . (151) 



The potential W{q) has the meaning of the free energy with the constraint that the overlap of our 
configuration r with the generic configuration a is equal to q. As far as we are interested in studying 
the g-dependance of W{q), we can set conventionally W{0) = 0. 

If one computes the functions r(e) and W{q) in a mean field model, one typically finds that 
the shape of the function W is characteristic of a mean-field system undergoing a first order phase 
transition. At high enough temperature W is an increasing and convex function of q with a single 
minimum for q = 0. Decreasing the temperature below a value Tf, where for the first time a point qf 
with W"{qf) = appears, the potential looses the convexity property and a phase transition can be 
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Figure 9: Phase diagram in the T — e plane. At the upper curve the low q solution disappear, at the 
lower curve the high q solution disappear and two locally stable solutions are present only in the region 
among the upper and lower curves. The middle curve the coexistence line where the two solutions 
have equal free energy. The coexistence line touches the axes e = at T = Tg, while the lower curve 
touches it at T = Tq. 

induced by a field. A secondary minimum develops at T^, the temperature of dynamical transition [7], 
signaling the presence of long-life metastable states. The height of the secondary minimum reaches 
the one of the primary minimum at T = Tg and thermodynamic coexistence at e = takes place. This 
is the usual static transition. In figure 8 we show the shape of the potential in the various regions. 

Therefore the potential W{q) has usually a minimum at q = 0, where M^(0) = 0. It may have a 
secondary minimum at q = qo- We have a few different situations: 

• At T > Td the potential W{q) has only the minimum at q = 0. The quantity qo cannot 

be defined and no valley with the equilibrium energy are present. This is more or less the 

definition of the dynamical transition temperature Tp. A more careful analysis [42] shows that 

for To < T < Ty there are still valleys with energy less than the equilibrium one, but these 

valleys cover a so small region of phase space that they are not relevant for equilibrium physics 
17 

• Exactly at T = Tq we sit at a phase transition point where some susceptibilities are divergent. 
This fact implies (in short range models) that there is a divergent dynamical correlation length 
that is related to dynamical heterogeneities [76]. 

• W{qD) > 0. This happens in an intermediate temperature region, above Tc, but below To, 
where we can put one replica a at equilibrium and have the second replica r in a valley near it. 
It happens that the internal energy of both the a configuration (by construction) and of the r 
configuration are equal to the equilibrium one. However the number of valley is exponentially 
large so that the free energy a single valley will be higher that the total free energy. One finds 

^'^In the REM limit (p — > oo) the temperature To goes to infinity. In this hmit the region T > Td does not exist. 
Therefore the dynamical transition is a new feature that is not present in the REM. 
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Figure 10: The full line is the function W{q) computed in the mean field approximation. The dashed 
line is the correct result (Maxwell construction). 

in this way that W{qD) > is given by 

W{qn) = ^ - (152) 

where Me is the average number of the valleys having the equilibrium energy [40, 44] . 

• W{qD) = 0. This happens in the low temperature region, below Tc, where we can put two 
replicas both at overlap and at overlap qEA without paying any prize in free energy. In this 
case qo = qEA- 

For T < To there is a relation (eq. (146) ) among the entropy inside a valley the entropy of 
the systems and is the configurational entropy, or complexity, that is given by the value of W at the 
secondary minimum. This W contribution vanishes at Tc and becomes exactly equal to zero for T < Tc 
[7] . 

Although the behavior of this potential function is analogous to the one found in ordinary systems 
undergoing a first order phase transition the interpretation is here radically different. While in ordinary 
cases different minima represent qualitatively different thermodynamical states (e.g. gas and liquid), 
this is not the case here. In our problem the local minimum appears when ergodicity is broken, and 
the configuration space splits into an exponentially large number of components. The two minima are 
different manifestations of states with the same characteristics. The height of the secondary minimum, 
relative to the one at g = measures the free-energy loss to keep the system near one minimum of 
the free energy (in configurations space). This is just the complexity TS, i.e. the logarithm of the 
number of distinct valleys of the system. 

The equation dW{q)/dq = e may have two stable solutions (that correspond to a local minimum 
of W{q) — eq) only in the region of the T — e plane shown in fig. 9. At the upper and low curves one 
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Figure 11: The shapes of the function q{e) for T > T^. the full line is the correct result and the dashed 
line is the output of a mean field approximation. 

of the two solutions loose its stability and it disappears: these two curves are the equivalent of the 
spinodal lines in usual first order transition. The point where the lower curve crosses the axis e = is 
the dynamical transition [39]: only at lower temperatures the two systems may remain with an high 
value of the overlap without having a force that keeps them together (i.e. e = 0). On the contrary the 
static transition is characterized by the fact that the coexistence line touches the axis e = 0. 

General arguments tell us that the free energy is a convex function of the q, so that we the correct 
shape of the function W can be obtained by the Maxwell construction (see fig. 10). In order to se the 
consequences of this fact on the definition of the complexity we can try to consider the function q{e) 
for temperatures less than Td shown in fig 11. 

As can be seen from the figures, the point where we evaluate the complexity (i.e. e = and high 
q) is always in the metastable region for T > Tg where we equilibrium complexity is non-zero. This 
causes an intrinsic ambiguity in the definition of complexity because the free energy in not defined 
with infinite precision in the metastable phase. However we can use the fact that the free energy is as 
function of e near the discontinuity point to extrapolate the high e free energy in the metastable 
region. This ambiguity becomes smaller and smaller more we approach the static temperature (the 
amount of the extrapolation becomes smaller and smaller) and in general it is rather small unless 
we are very near to the dynamic phase transition. This ambiguity is not important from practical 
purposes; however it implies that there is no sharp, infinitely precise definition of the equilibrium 
complexity. If we forget this intrinsic ambiguity in the definition of the complexity we may arrive to 
contradictory results. 
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4.3 Complexity and replicas 

As we have seen we can write 

Z{(3) = J2 exp(- = f cW{f, (3) exp(-/3iV/) , (153) 

a •' 

where /a(/3) is the free energy density of the valley labeled by a at the temperature (3"^, and M{f,P) 
is the number of valleys with free energy density less than /. 

We have also seen that M(f, (3) = exp(A/'S(/, 13)), where the configurational entropy, or complexity, 
'E{f,/3) is positive in the region / > fo{(3) and vanishes at / = /o(/?)- The quantity fo{P) is the 
minimum value of the free energy: J\f{f,(3) is zero for / < fo{(3) [39, 44, 43]. 

If the equation 

has a solution at / = f*{(3) (obviously this may happens only for /*(/?) > /o (/?)), we stay in the liquid 
(high temperature) phase. Here the free density is given by 

feg = r-r'm*,p) (155) 

andE(r,/3) = S*(/3). 

Otherwise we stay in the glass (low temperature) phase and 

feg = MP) ■ (156) 

In order to compute the properties in the glass phase wc need to know (3): a simple strategy 
to compute the complexity is the following. We introduce the modified partition function 

Z{r, P) = eM-NjGir, P)) = J2 eM-lNfaiP)). (157) 

a 

It is evident that Z{P; (3) is the usual partition function and G{P\ 0) is the usual free energy. Using 
standard thermodynamical arguments it can be easily proven that in the limit N ^ oo one has: 

lG{r,P)=7f-mf), f = ^h^M)l. (158) 

The complexity is obtained from G{'y; P) in the same way as the entropy is obtained from the usual 
free energy [40, 43]: 

E(A/) = 5^. (159) 

A few observations are in order: 

• In the new formalism 7, the free energy and the complexity play respectively the same role of 
(3, the internal energy and the entropy in the usual formalism. 

• In the new formalism f3 only indicates the value of the temperature that is used to compute the 
free energy and 7 controls which part of the free energy landscape is sampled. 
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• When j3 —> CO (at least in mean field models) we sample the energy landscape: 

Z(7;oo) = ^exp(-7A/"ea) = / z/(e)deexp(-7A/"e) (160) 

a ■' 

where Ca are the minima of the Hamiltonian and v{e) the density of the minima of the Hamil- 
tonian. 

• The equilibrium complexity is given by S*(/3) = S(/3;/3). On the other hand S(7;oo) give us 
information on the minima of the Hamiltonian. 

In principle it is possible to get the function S(/) by computing directly the number of solution 
of the TAP equations for a given value of the free energy density. However it is simpler to obtain it 
by using the replica formalism and it is reassuring that one gets the same results with both methods 
[40, 39, 44, 45, 46, 47]. 

The computation of the modified partition function Z{'-^; (3) can be easily done in the replica 
formalism [39, 44]. If we consider a system with m replicas (with m integer) and we constrain them 
to stay in the same state we find that 

Z{(3, m) = Y, eM-PmNfaiP)) (161) 

a 

This expression coincide with Z^^y; (3) for 7 = m(5. Therefore there is a very simple way for computing 
G'(7;/?)). We must consider the partition function of m replicas that are constrained to stay in the 
same state, i.e. there are at a large value of q where q is chosen in a self consistent way.. 

Let us firstly see how this approach works in the REM. In the REM the configurations and the 
states practically coincide: there is no (5 dependence of the complexity. The REM is defined by 
the property that S(e) = — ln(2), so that the computation of S(e) using the replicas may look 
pointless, however it is instructive to illustrate the point. 

The contribution to the partition function of m replicas coming from the region of phase space 
where all the m configuration are identical is given by 

Z(/3, m) = 2^ exp(- ^ {Pmf) (162) 

We thus get 

7G(7;/?)=ln(2)-i72 (163) 

and from the previous equation we can read back the expression for the complexity. 

In the p-spin model we have to find out the partition function of a systems where all the m replicas 
are in the same block and are characterized by an high value of q. Here, given the value of m, we have 
to look for a solution q* of the equation for q of the replica approach: 

dF 

^=0, (164) 

with q* 7^ 0, where the potential F{q, m) is the replica potential introduced in the previous section. 
We thus find that 

G{(3m-(3))=F{m,q*) (165) 
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The computation that we have done before for the statics contains the whole information needed to 
compute also the complexity and some other properties of the metastable states. Indeed we get 



= ^^^^^ . (166) 
m 

and f{m) can be obtained by Legendre transform or by using the relation 

f{m) = g^-?'^*^ = Q*) - mE{m) . (167) 

One finally finds the complexity as function of the free energy, by eliminating m. 

S(/) = S(m(/)) (168) 

The dynamical temperature is is the highest temperature where eq. eq. (164) has a solution at 
near m = 1, while the static critical temperature is the first temperature where 

F{q{m),m)\m=o = F {q{m) , m)\m=i (169) 

The equilibrium transition temperature is given by the condition that the equilibrium complexity 
satisfies the condition 

E..5q!^U., = o. (170) 

In this way we have recovered the results of one step replica symmetry breaking (together with 
the mysterious condition dF{m)/dm = 0) from general principles. 

Although we have based our discussion on mean-field model, we expect that the qualitative features 
of the phase diagrams presented survive in finite dimension. We believe that the existence of a 
coexistence line, terminating in a critical point, is a constitutive feature of systems whose physics 
is dominated by the existence of long lived metastable states like glasses. These predictions of can 
be submitted to numerical test in glassy model systems as like e.g. Lennard-Jones or hard spheres, 
or polymer glasses. For example the identification of the complexity S as the free energy difference 
between the stable and the metastable phases allows an other way to measure of this quantity in a 
simulation. Indeed the ending of the transition lines in a critical point implies that the metastable state 
can be reached via closed paths in phase diagram leaving always the system in (stable or metastable) 
equilibrium; the free energy difference of the two phases can be computed integrating the derivative 
of the free energy along such a closed path. 

We have to study the shape of the function F(to, q). At fixed m as function of q it may have one 
of the forms shown in fig. 8. Below the dynamical transition near m = 1 it has the shape of the 
lower curve fig. 8. By decreasing m the shape of this function modifies and the secondary minimum 
disappears. There is a temperature-dependent region of m where the equation dF/dq = has a 
solution at non-zero q: in this region we can compute the complexity. It is possible that if we compute 
the small fluctuations in this region using the techniques of the next section we finds for some values 
of m a not consistent result, i.e. e negative spectrum. This phenomenon may indicate that the allowed 
m region is smaller than that indicated by the condition of the existence a solution to the equation 
dF/dq = 0. However it is also possible that mor complex phenomena are present, that are not fully 
understood at the present moment. 
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4.4 A summary of the results 

I will now summarize the results. As we have seen we can distinguish a few temperature regions. 

• For T > Tf the only minimum of the free energy functional is given by high temperature result: 
we call it the liquid minimum (in the spins models described above it has to zero magnetization). 

• For Tf > T > To there is an exponentially large number of minima [23, 42, 44]. For some values 
of the free energy density the complexity S is different from zero, however the contribution to 
the free energy coming from these minima is higher that the one coming from the liquid solution 
with zero magnetization. As discussed also in Cugliandolo's lectures the value To coincides with 
the critical temperature of the mode coupling approach and in the glass community is called Tc. 
The real critical temperature of the model, that we have called Tc up to now is called Tk for 
reasons that will be clear in the next section. 

• The most interesting situation happens in the region where Td > T > Tc {or Tc > T > Tk 
using the glassy notation). In this region the free energy is still given the high temperature 
solution (with zero magnetization in spin models), It is extremely surprising [38, 40] that the 
free energy can be written also as the sum of the contribution of an exponentially large number 
of non-trivial minima as in eq. 142).. 

Although the free energy is analytic at To , below this temperature the the system at each given 
moment may stay in one of the exponentially large number of minima. The time (r) to jump 
from one minimum to an other minimum is quite large and it is controlled by the height of the 
barriers that separate the different minima. In the mean field approximation (i.e. for infinite 
range models) it is proportional to exp{AN) with non-zero A. In short range models at finite 
dimensions we expect that the barriers are finite and r ~ tq exp(/3A(T)). The quantity /3A(T) 
is often a large number also at the dynamical temperature [48] (e.g. 0(10)) and the correlation 
time will become very large below To and for this region To is called the dynamical transition 
point. The correlation time (that should be proportional to the viscosity) should diverge at Tk- 
The precise form of the this divergence is not well understood. It is natural to suppose that we 
should get divergence of the form exp(y4/(T — TkY) for an appropriate value of v [49], whose 
reliable analytic computation is lacking [7, 37]. The value v = 1 (i.e. the Vogel Fulcher law) is 
suggested by the experiments. 

The equilibrium complexity is different from zero (and it is a number of order 1) when the 

temperature is equal to T^ and it decreases when the temperature decreases and it vanishes 
linearly at T = Tk- At this temperature (the so called Kauzmann temperature) the entropy of 
a single minimum becomes equal to the total entropy and the contribution of the complexity to 
the total entropy vanishes. At an intermediate temperature Tg the correlation time becomes so 
large that it cannot be observed any more by humans. 

• In the region where T < Tk the free energy is dominated by the contribution of a few minima 

of the free energy having the lowest possible value. Here the free energy is no more the analytic 
continuation of the free energy in the fluid phase. A phase transition is present at Tk and the 
specific heat is discontinuous here. 

4.5 Some consideration on the free energy landscape and on the dynamics 

The free energy landscape is rather unusual; we present the following pictorial interpretation (fig. 12), 
that is a rough simplification [50] . At a temperature near to the system stays in a region of phase 
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Figure 12: The qualitative dependence of the free energy as function of the configuration space in the 
region relevant for the dynamical transition, i.e. for T < Tjj. 

space that is quite flat and correspond of a minimum of the total free energy. On the contrary below 
To the phase space is similar to the one shown pictorial in fig. 12. The region of maxima and minima 
is separated by the region without barriers by a large nearly flat region. The minima in the region at 
the left are still present also when Tf > T > To, but they do not correspond to a global minimum. 

At temperatures higher than To the system at thermal equilibrium stays in the right portion of 
fig. 12. When the temperature reaches To the system arrives in the flat region. The flatness of the 
potential causes a Van Hove critical slowing down that is well described by mode coupling theory 
[12, 18] (that is exact in the mean field approximation). 

In the mean field approximation the height of the barriers separating the different minima is 
infinite and the temperature To is sharply defined as the point where the correlation time diverge. 
The precise meaning [39] of the dynamical temperature beyond mean field approximation has already 
been discussed. 

Let us start from a very large system (of N particles) at high temperature and let us gradually 
cool it. It would like to go at equilibrium in the region with many minima. However coming from high 
free energy (from the right) it cannot enter in the region where are many maxima; if we wait a finite 
amount of time (the time to crosses the barriers diverges as exp(^A^). the system remains confined 
in the flat region. In this case [9, 10] the so called dynamical energy, 

Ed = lim lim E{t,N), (171) 

is higher that the equilibrium free energy. The situation is described in fig. 4. 5. The difference of the 
static and dynamic energy is an artifact of the mean field approximation if we take literarily the limit 
t — > oo However it correctly describe the situation on laboratory times, where metastable states are 
observed. 

In the mean field approximation very interesting phenomena happen below when the system 
is cooled from the high temperature phase due to the fact that the system does not really go to an 
equilibrium configuration but wanders in the phase space never reaching equilibrium. The phenomena 
are the following: 
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Figure 13: The qualitative behaviour of the equihbrium energy and of the dynamical energy as function 
of the temperature. 

• The energy approaches equilibrium slowly when the system is cooled from an high energy 
configuration [9, 51]: 

E{t, T) = Ed{T) + B{T)t~^^^\ (172) 

where the exponent A(T) does not vanish linearly at zero temperature as happens for an activated 
process. 

• Aging is present, i.e. the correlation functions and the response functions in the region of large 
time do depend on the story of the system [17, 16, 10]. 

• In the region where aging is present the fluctuation dissipation theorem is no more valid. New 
generalized relations are satisfied [9, 10, 52, 53, 54, 55], that replace the equilibrium fluctuation 
dissipation theorem. 

These phenomena will be discussed in details in Cugliandolo's lectures. 

It is not clear how to compute in general the quantity E£,[T). There is a very simple recipe: Ejy is 
the largest energy where > 0. According to that recipe one has to look to the smallest value of 

m where the equation for q has a consistent solution (as discussed in the previous section). This last 
condition is equivalent to impose that the spectrum of small fluctuations of the free energy (defined in 
the next section) has a gap. The value of E£) is characterized by the fact the replica broken solution 
becomes unstable aX E > Eq (marginal stability). This marginalistic approach correctly gives Eq in 
the p-spin spherical model, it is not clear if it true in general [42]. 

4.6 Small fluctuations 

We have see that in the real world glassy systems have only one transition with divergent correlation 
time (at temperature Tr-). However in the idealized world of mean field theories there is a second 
purely dynamics transition To at higher temperatures [7]. As it happens in many cases, slow relaxation 
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is related to the existence of zero energy modes and this statement is true also here. This statement 
can be easily verified in spin models where the mode coupling theory is exact and simple computations 
are possible. 

In spin models we concentrate our attention on the Hessian of the free energy, defined as 

M(.,fc) = ^'.y , (173) 
om{i)om[k) 

where F[m\ is the (TAP) free energy as function of the magnetization and the magnetizations do 
satisfy the stationarity (TAP) equations: 

Performing the appropriate computations [56, 57] we finds that the spectral density of the Hessian in 
these infinite range models has always a semicircular form: 



p(A) oc ^{\-y^{T)){\-v{T)) , (175) 

typical of random matrices [58]. 

• At temperatures T > r£))there are no non-trivial thermodynamically relevant solutions of the 

equation eq. (174) , however the dynamics is dominated by quasi-solutions of the previous 
equations, i. e. by magnetizations such that the left hand side of the previous equation is not 
zero, but small [59].. The Hessian M of the quasi-solutions has negative eigenvalues and its 
spectrum has qualitatively the shape shown in fig. (14). These quasi stationary points of F look 
like saddles. 

• At the dynamical transition point T = Td the quasi stationary points becomes real solutions of 

the equations (174). They are essentially minima: the spectrum of the Hessian is non-negative 
and it arrives up to zero. As it can be checked directly, the existence of these zero modes is 
responsible of the slowing down of the dynamics. The different minima are connected by flat 
regions so that the system may travel from one minimum to an other [50]. 

• At low temperature the mimima become more deep, the spectrum develops a gap as shown in 
fig. 14 and the minima are no more connected by flat regions. In mean field models the system 
would remains forever in one of these minima. If the system starts from an high temperature 
configuration it cannot reach these configurations. 

This picture is not so intuitive because it involves the presence of saddles with many directions in 
which the curvature is negative, and it is practically impossible to visualize it by making a drawing 
in a two or a three dimensional space. 

This qualitative description can be easily verified in models where the mean field approximation 
is exact. However, if we try to test it in finite dimensional models, we face the difficulty that the 
free energy functional F[m] is a mythological object whose exact form is not known and consequently 
the eigenvalues of its Hessian cannot be computed. A n alternative approach consists in studying the 
properties of the so called instantaneous normal modes (INM) [60, 57, 61] and by the saddle normal 
modes (SNM) [62, 63, 64, 48]. For reasons of space this interesting point cannot be discussed any 
more the reader is invited to look to the original literatures. 
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Figure 14: The qualitative behaviour of the spectrum in mean field approximation above Tc, (full 
line), at Tc (dot-dashed line) and below Tc (dashed line) as function of the eigenvalue A). 

These computations are particular relevant for glasses where the low temperature spectral density 
can be experimentally measured and one finds the famous Boson peak In order to explain the 
Boson peak one has to do realistic computations, where one has to take into account the spectrum 
of phonons and this has been done or using the mode coupling formalism [67], or a microscopical 
approach [68]. 

5 Structural relations 

In this section we shall define a new equilibrium order parameter function P{q), that is connected 
also to fluctuation-dissipation ratio X{q) (see Cugliandolo's lectures). This can be done by studying 
of the linear response to some special sets of perturbations of the original Hamiltonian [69]. This 
method has been recently used to derive interesting properties of the overlap distribution at equilibrium 
[70, 71, 72, 73]. 

5.1 Stochastic stability 

Our aim it to prove that the moments of P{q) are related to the response of the system when one adds 
an appropriate perturbation [70]. This approach allows us to define all the relevant quantities in the 
case of single large system (in the infinite volume limit), while in the previous approach the function 
P{q) was defined as the probability distribution in an ensemble of different systems, characterized by 

^^The Boson peak is defined as a bump at some small but non-zero value of the eigenvalue z of the spectral density p{z) 
divided by the Debye density of states z'^) that has been observed in many material [65] and in numerical simulations 
[66]. 
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different realizations of the disorder. Tliis difference is crucial if we consider the case (like glasses) 
where no disorder is present 

In the case of spin systems an appropriate form of the long-range perturbation s given by: 

1,N 

Hp{(^) = Kh,...,ip(^ii ■■■(^ip, (176) 

ii,ip=l,N 

where the couplings Ki^^,,,^i^ are independent Gaussian variables with zero mean and variance Kf^ ^ = 
l/(27VP-i). 

The Hamiltonian is 

H, = Hi + eHp, (177) 
The canonical average of Hp verifies, for all values of e, the relation 

JH^ = -PeN (l-Jdq PM q^^ , (178) 

irrespective of the specific form of Hi. Here the function Pe{q) is the probability distribution of the 
overlap q in the presence of the perturbing term in the Hamiltonian; the average is done over the new 
couplings J at fixed Hj. The derivation, involves only an an integration by parts in a finite system. 

However the previous equation is strange. The function P^{q)\e=Q depends on the instance of the 
problem also in the infinite volume limit, while, for e 7^ 0, {H^^) is a thermodynamic quantity shat 
cannot fluctuate in the infinite volume limit when we change the instance of the system, (at least 
for generic e). A further difficulty appears when there is an additional symmetry in absence of the 
random perturbation that is broken by the random perturbations (e.g. the spin reversal symmetry 
for a spin system in absence of a magnetic field or the translational invariance for glasses). Therefore 
for a large system we may expect that 

P{q) / Pi{q), (179) 

where Pi{q) where is the usual overlap probability distribution computed at e = 0, that may depend 
on the instance, while \P{q) is the limit e — of the function Pe{q) (where the limit is computed 
outside the cross-over region, i.e. e >> N^^"^) and it does not depend from the instance. 

We need to understand better how the equilibrium state in the presence of the perturbation eHp is 
related to the equilibrium state at e = 0. Some complications arise whenever the equilibrium expecta- 
tion value of Hp is not the same for all pure phases of the unperturbed systems. Then the limit value of 
iip as e — will be the one corresponding to the favored phases. A simple example is the Ising model 
in the ferromagnetic phase, where one adds a negative magnetic field term as a perturbation. Also 
in the vanishing field limit, the the system stays in a state with negative spontaneous magnetization, 
while the unperturbed measure corresponds to a mixture of the positive and negative magnetization 
pure states. 

In presence of many equilibrium states that enter in the Gibbs-Boltzmann measure, as it happens 
where replica symmetry is broken, the situation is rather complex. Indeed, the stochastic perturbations 
that we have considered will in general reshuffle the weights of the different ergodic components in 
the Gibbs measure, or even change their nature, and this changes the P{q) function to a different one 

(^(g))- 

^^In a glass we still have the possibility of averaging over the total number of particles. 
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The principle of stochastic stabiUty assumes that if consider an appropriate ensemble as for the 
initial random system we have that 

P{q) = = Hq) (180) 



where here the bar denotes the average over the instances of the system in the appropriate ensemble. 

An intuitive motivation for assuming stochastic stability is that perturbations wc use are random 
and they are not correlated with the original Hamiltonian. So they should change the free energies of 
the various pure states of the original systems by random amounts. Stochastic stability assumes in 
glassy systems the distribution of these free energies is stable under independent random increments, 
as has been shown in mean field (in fact this property lies at the heart of the cavity method [2]). If 
this is the case then the two functions P{q) and P{q) will coincide. We could also say that the equality 
of P{q) and P{q) indicates that the systems responds to a random perturbation as a generic random 
system. 

There are cases where stochastic stability trivially fails, i.e. when the original Hamiltonian has an 
exact symmetry, that is lifted by the perturbation. The simplest case is that of a spin glass with a 
Hamiltonian invariant under spin inversion. In this case P{q) = P{—q), since each pure state appears 
with the same weight as its opposite in the unperturbed Gibbs measure. On the other hand, if we 
consider Hp with odd p, this symmetry is lifted. This means that in the e ^ limit only half of 
the states are kept. If the reshuffling of their free energies is indeed random, then we shall have 
P{q) = 20{q)P{q) = P{q). The same type of reasoning applies whenever the overlap q transforms 
according to a representation of the symmetry group of the unperturbed Hamiltonian Hq. 

If this trivial effect of exact symmetries is taken into account, for a large class of systems, the 
function P{q) in the limit of small perturbations tends to the order parameter function P{q) of the 
pure system where the exact symmetries are lifted. This continuity property is called stochastic 
stability. Ordinary systems without symmetry breaking and mean-field spin glasses are examples 
of stochastically stable systems. In ergodic systems, the equality of P and P is immediate, both 
functions consist in a single delta function. Thus, the problem of deriving the equality between P 
and P, appears only when there are coexisting phases are unrelated by symmetry. Unfortunately, 
we are not able to characterize the class of stochastically stable systems. In particular there is no 
rigorous proof that short-range spin glass, for which our theorem is most interesting, belong to this 
class. However, stochastic stability has been established rigorously in mean field problems [70, 71]. 

If one studies more carefully the problems, one finds that stochastic stability has far reaching 
consequences, e.g. 

Pl{qi)Pl{q2) = lPiqi)P{q2) + lp{qi)S{qi - qi) (181) 

These (and other) relations have been carefully numerically verified in also in numerical simulations 
of three dimensional spin glasses models [53] strongly suggesting the validity of stochastic stability. 

5.2 A simple consequence of stochastic stability 

Stochastic stability is a very powerful property, and it is the ingredient that allows to relate the prop- 
erties of the low lying configurations, that dominate the Gibbs measure, to those of the configurations 
much higher in energy that are seen in the dynamics. This is most easily explained in the usual 
framework of replica-symmetry breaking, considering an approximation with only two possible values 
of the overlap, qo among different states and qi among the same state (i.e., one-step replica-symmetry 
breaking). The probability of finding a state with total free energy = F is given by p{AF), where 
AF = F — Fq and Fq is a reference free energy the equilibrium free energy. The weight of each state 
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is given by 

wa oc exp(-/3Fa). (182) 

In one-step replica-symmetry breaking, the states that contribute to the Gibbs measure have nearly 
degenerate free energies. We have already seen the non-extensive fluctuation of their free energies, 
corresponding to the low AF regime of p{AF), is given by [25, 92] 

p{AF) oc exp(/3mAF), (183) 

and the function P{q) is given by 

P{q) = mS{q - qo) + (1 - m)Siq - qi) . (184) 

Stochastic stability forces the function p{AF) to be of the form (183), not only when AF is finite, 
but also in the range where AF is extensive but small (say of order eA'^). Indeed it imposes that the 
form of the function p{AF) remains unchanged (apart from a possible shift in Fq) when one adds a 
small random perturbation [69]. 

Let us consider the effect of a perturbation of strength e on the free energy of a state, say a. The 
unperturbed value of the free energy is denoted by F^. The new value of the free energy Gq, is given 
by Ga = Fa + era where are identically distributed uncorrelated random numbers. Stochastic 
stability implies that the distribution p{G) is the same as p{F). Expanding to second order in e we 
see that this implies dp/dF oc d^p/dF'^, whose only physical solution (apart the trivial one p{F) = 0, 
that corresponds to non-glassy systems) is given by eq. (183) We see that stochastic stability fixes 
the form of the function p and therefore connects in an inextricable way the low and the high free 
energy part of the function p. 

This remark explain how it is possible that stochastic stability tells us something on the dynamics 
in the aging regime (as we shall see in the next section). In the dynamical evolution from an higher 
temperature initial state, the difference between the total free energy at time t and the equilibrium 
value will be always of order A'', with a prefactor going to zero when t goes to infinity. One could argue 
that the d3mamics probes the behavior of the function p{AF) at very large argument, and should not 
be related to the static property that depend on the function p for small values of the argument. 
However stochastic stability forces the function p{AF) to be of the form eq. (183) , also in the range 
where AF is extensive but small, and this objection is no more valid. 

5.3 Fluctuation dissipation relations 

Let us now discuss the case of dynamics of a systems that at times starts from a non-equilibrium 
configuration. Here the relevant quantities are the two times correlations functions and the response 
functions, e.g. the correlation defined as 

Cft^.t)^ ^'-''"";""'"'"' (185) 

The finite-time response and correlation functions involved in the definition of the FDR are con- 
tinuous functions of e for e ^ 0. It is not evident if the limit is uniform in time, i.e. if the infinite 
times and the e ^ limits do commute. The linear response regime may shrinks to zero as the time 

^''The same conclusion could be obtained using the methods of reference [73] computing the sample-to-sample fluctu- 
ations of the function Pj{q), that in this case, where ultrametricity is trivially satisfied, are completely determined by 
the knowledge of of the function P{q). 
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goes to infinity. This possibility shows up when the perturbations favor one phase (as discussed in the 
introduction). However here we consider a random perturbations and the expectations of Hp vanishes 
at e = 0: it is reasonable to assume that the linear response regime survives at very long times 

If the infinite times limit and the e ^ limit do commute (i.e. the dynamical form of the stochastic 
stability), the dynamics is strongly constrained. For example, if we consider the function X{q) that 
parametrizes the violations of the fluctuation dissipation relations in off-equilibrium experiments [9, 
10], the previous assumptions imply that: 

= f . (186) 

. The proof is simple it only involves integrations by parts. The details can be found in the original 
papers [69]. 



6 A short introduction to glasses 

Glasses [74] are roughly speaking liquids that do not crystallize also at very low temperature (to be 
more precise: glasses also do not quasi-crystallize). 

These liquids can avoid crystalization mainly for two reasons: 

• The liquid docs not crystallize because it is cooled very fast: the crystallisation time may become 
very large at low temperature (e.g. hard spheres at high pressure). The system should be cooled 
very fast at temperatures near the melting point; however if crystalization is avoided, and the 
temperature is low enough, (e.g. near the glass transition) the system may be cooled very slowly 
without producing crystalization. 

• The liquid does not crystallize even at equilibrium. An example is a binary mixture of hard 

spheres with different radius: 50% type A (radius r^, 50% type B (radius rA, where R denotes 
tb/ta- If -77 < i? < .89 (the bounds may be not precise), the amorphous packing is more dense 
than a periodic packing, distorted by defects. 

Which of the two mechanism is present is irrelevant for understanding the liquid glass transition. 
Other examples of glassy systems are binary mixtures: in this case we have 2 kinds of particles 
and the Hamiltonian of N particles is given by 

^= E E E ya,b{Xa{i)-X,{k)) , (187) 
a,6=l,2 i=l,Ar(a) k=l,N(b) 

where a = 1, 2; N{a) = Nc{a); J2a=i,2 c(«) = 1- In this case the values of the N concentrations c and 
the 3 functions Va,b{x) describe the model. Well studied case are: 

• A power potential (non-realistic, but simple), e.g. c(l) = .5, c(2) = .5, Va^b(x) = RafiX~^'^. 

• Lennard-Jones potential (more realistic), e.g. c(l) = .8, c(2) = .2, Va^b{x) = Ra,bX^^^ — AafiX^^. 

There are some choices of the parameters that are well studied such that the system does not 
crystallize. One has been introduced by Kob and Anderson in the L-J case [75]: it corresponds to a 
particular choice of the parameters R and A. 

^^Let us stress that the existence of a hnear response regime uniform in time a question susceptible of experimental 
investigation. 
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Figure 15: The viscosity as function of the temperature according to the Vogel Fulcher law (Tx = 1 
full line) and to the mode coupling theory [Tk = 1.375 dotted line). 

There are many other material that are glass forming, e.g. short polymers, asymmetric molecules 
(e.g. OTP) . . . 

The behaviour of the viscosity in glass forming liquid is very interesting There are two regimes: 

• In an high temperature region the mode coupling theory [12] is valid: it predicts r] oc (T — T(i)~^ , 
where 7 is not an universal quantity and it is 0(1). 

• In the low temperature region by the Vogel Fulcher law [49] is satisfied: it predicts that r] oc 
exp{A{T — Tk)~^)- Nearly tautologically fragile glasses can be defined as those glasses that have 
Tk 7^ 0; strong glasses have Tk ~ 0. 

At the glass temperature (Tg) i the viscosity becomes so large that it cannot be any more measured. 
This happens after an increase of about 18 order of magnitude (that correspond to a microscopic time 
changing from 10~^^ to 10^ seconds): the relaxation time becomes larger than the experimental tine. 

A characteristic of glasses is the dependance of the specific heat on the cooling rate. There is a 
(slightly rounded) discontinuity in the specific heat that it is shifted at lower temperatures when we 
increase the cooling time. 

For systems that do crystallize if cooled too slowly, one can plot AS = S (liquid) — S (crystal) 
versus T in the thermalized region. One gets a very smooth curve whose extrapolation becomes 
negative at a finite temperature. A negative AS" does not make sense, so there is a wide spread belief 
that a phase transition is present before (and quite likely near) the point where the entropy becomes 
negative. Such a thermodynamic transition (suggested by Kaufmann) would be characterized by a 

^■^The viscosity can be defined microscopically considering a system in a box of large volume V: T^^u(t) is the total stress 
tensor at time t. We define the correlation function of the stress tensor at different times: (r^,^(t)rp,o-(0)) = V^5'^,i/,p,(T(i) 
Neglecting indices, 11 oc J dtS{t) ~ t~", where r is the characteristic time of the system (in a first approximation we can 
suppose that a = 1). 
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Figure 16: The specific heat and the entropy excess AS* for various coohng rates. 

jump of the specific heat. Quite often this temperature is very similar to the temperature found by 
fitting the viscosity by the Volker Fulcher law and the two temperature are believed to coincide. 

Usually a second order a thermodynamic phase transition induces a divergent correlation time. 
However the exponential dependance of the correlation time on the temperature is not so common (in 
conventional critical slowing down we should have a power like behaviour); an other strange property 
of the glass transition is the apparent absence of an equilibrium correlation length or susceptibility 
(linear or non- linear) that diverges when we approach the transition point^'^. 



7 The replica approach to structural glasses: general formalism 

In this section we write down the formulas corresponding to the replica approach introduced in the 
previous section. We keep here to the case of simple glass formers consisting of N particles interacting 
by a pair potential v{r) in a space of dimension d. 

The reader may notice that in all the example that we have considered up to now a quenched 
disorder was present; this feature is not present in glasses, where no random quenched variables are 
present in the Hamiltonian. However this in not a difficulty as far as there is no need of a quenched 
random disorder to use the replica formalism as it is clear from the previous sections. However at 
the beginning it was believed that a random quenched disorder was necessary in order to use the 
replica approach. Only much later [13, 77, 78, 79] it was realized that the replica method could be 
used in translational invariant models where no disorder is present; some of these models behaves in 
a way very similar to real glasses: one can define both To and Tk and at low temperature the system 
crystallizes. 

■^•^There is a dynamical correlation length that diverge at the critical temperature that is observed numerically and it 
is predicted theoretically [76]. 
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7.1 The partition function 

The usual partition function, used e.g. in the hquid phase, is in tree dimension 

1 r ^ 



=1 



We wish to study the transition to the glass phase through the onset of an off-diagonal correlation in 
replica space[80, 81]. We use m replicas and introduce the Hamiltonian of the replicated system: 

m 

Hm= E E^«-4) + ^E E ^«-^') (189) 

l<i<j<Na=l i=lNa<b=m 

where w is an attractive interaction. The precise form of w is unimportant: it should be a short range 
attraction respecting the replica permutation symmetry, and its strength e that will be sent to zero in 
the end. For instance one could take ^ 

= (^) ^^^^^ 

with c is of the order of 0.2 times the typical intermolecular distance. A positive value of e forces all 
the particles one near to the other ones. If two systems stay in the same state, the expectation value 

of w represents the self overlap and it very near to 1. 
The partition function of the replicated system is 



1 „ N rn 

Zm^j^ \[\[{d'xt)e-P"- (191) 

1=1 a=l 



7.2 Molecular bound states 



At low enough temperature, we expect that the particles in the different replicas may stay close to 
each other. The role of the attractive term w is to insure that all replicas fall into the same glass 
state,: also for small e the particles in different replicas stay at the same place, apart from some 
thermal fluctuations: 

Thermal fluctuations are relatively small throughout the solid phase (one can see this from the 
Lindeman criterion) and diffusion is very small, one can identify the molecules and relabel all the 
particles in the various replicas in such a way that the particle j in replica a always stays close to 
particle j in replica h. All the other relabelings are equivalent to this one, producing a global factor 
jYjm-i ^j-^ ^YiQ partition function. 

We therefore need to study a system of molecules, each of them consisting of m atoms (one atom 
from each replica). It is natural to write the partition function in terms of the variables that 
describe the centers of masses of the molecules, and the relative coordinates u", with = Ti + u^ and 
E„< = 0: 

N N m N 



= /^/n(^Snn('^v)nh'^(E< 



i=l i=l a=l i=l 



exp I -p E ^i^i - +<-«")- E E Wi^A) J (192) 

i<j,a i a,b 
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7.3 The small cage expansion 

In order to transform these ideas into a tool for doing explicit computations of the thermodynamic 
properties of a glass [81] we have to use an explicit method for computing the free energy as function 
of the temperature and m. As is usually the case, in the liquid phase exact analytic computations are 
not possible and we have to do some approximations. In this section wc shall use the fact that the 
thermal fluctuations of the particles in the glass are small at low enough temperature: the size of the 
'cage' seen by each particle is therefore small, allowing for a systematic expansion. What we will be 
describing here are the thermal fluctuations around the minimum of the potential of each particle, in 
the spirit of the Einstein model for vibrations of a crystal. 

We start from the replicated partition function Zm described in molecular coordinates in (192). 
Assuming that the relative coordinates are small, we can expand w to leading order and write: 

, „ TV N m N / \ 

■ '' i=l 1=1 a=l i=l\ a J 

exp I -/3 - + < - ^j) - ^ E E« - ^'n (193) 

y i<j,a i a,b J 

In the end we are interested in the limit e = (1/a) — > 0. We would like first to define the size A of the 
molecular bound state, that is also a measure of the size of the cage seen by each atom in the glass, 

^ ' ' i a,b 

{d is the dimension that we have taken equal to 3 and N is the number of particles). We Legendre 
transform the free energy (j){m,a) = —{T/m) log Zm, introducing the thermodynamic potential per 
particle t/j^m, A): 

^{m, A) = (t){m, a) + Td^^-^ - (195) 

2 a 

What we want to see is whether there exists a minimum of at a finite value of A. 

At low temperatures, this minimum should be at small A, and so we shall seek an expansion of 
'0 in powers of A. It turns out that it can be found by an expansion of (j) in powers of a, used as an 
intermediate bookkeeping in order to generate the low temperature expansion. 

This may look confusing since we are eventually going to send a to oo. However this method 
is nothing but a usual low temperature expansion in the presence of an infinitesimal breaking field. 
For instance if one wants to compute the low temperature expansion of the magnetization in a d- 
dimensional Ising model in an infinitesimal positive magnetic field h, the main point is that the 
magnetisation is close to one. One can organise the expansion by studying first the case of a large 
magnetic field, performing the expansion in powers of exp(— 2^), and in the end letting ^ — ^ 0. A 
little thought shows that the intermediate -large h- expansion is just a bookkeeping device to keep 
the leading terms in the low temperature expansion. What we do here is exactly similar, the role of 
h being played by 1/a. 
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7.3.1 Zeroth order term 
We use the equivalent form: 

In the Umit a — > 0, the identity 

m 



-, . N m fly I 

Zmioc) = ^ / n n (^'o n e ^« - ^.") - ^ E(^" - 

j=l a=l i 4/^^ \ i<i\o j.o 



(196) 



impUes that: 

Z^(a) = f--) _/jJdx,exp -^m^.;(X,-X,) . (198) 

In this expression we recognise the integral over the Xj's as the partition function ZuqiT*) of the 
liquid at the effective temperature T*, defined by 

T* = T/m . (199) 

Therefore the free energy, at this leading order, can be written as: 

(3<tP{m, a) = log ^ - A log(^) _ J_ bg Zu,{T*) (200) 

The result is intuitive: in the limits where the particles of different replicas stay at the same point, 
the Hamiltonian for m replicas is just the usual one, multiplied by m. 

7.3.2 First order term 

In order to expand to next order of the expansion, we start from the representation (193) and 
expand the interaction term to quadratic order in the relative coordinates: 

exp f -f E EK - - u^ f"%"'^ - ^ E« - y^f ' 

(where for simplicity we have not introduced the indices ^ and running from 1 to d, that denote 
space directions). Notice that in order to carry this step, we need to assume that the interaction 
potential v{r) is smooth enough, excluding hard cores. 

After some computations one finds that the free energy to first order is equal to: 

P4>{m, a) = 

3(m — 1) , 1 3(1 — m) , 27r 3 , 1 , ^ 
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where the constant C is proportional to the expectation value of the Laplacian of the potential, in the 
liquid phase at the temperature T*: 



C^W-^T.i^<'^-'3)Y (202) 
Differentiating the free energy with respect to l/a gives the equation for the size of the cage: 

Expanding this equation in perturbation theory in A we have: 

a = mA- -i- —A^ (204) 

3(m - 1) ^ ' 

The Legendre transform is then easily expanded to first order in A: 

X'^-m) A 



/3tp{m, A) = P(f){m, a) + 3- 



2 a 



= ^ log(27r^) - PmAC + !iL_Z!^ _ J_ log m - 4- log Zu,{T*) (205) 

2m 2m 2m mN 

This very simple expression gives the free energy as a function of the number of replicas, m, and 

the cage size A. We need to study it at m < 1, where we should maximise it with respect to A and m. 
The fact that we seek a maximum when m < 1 instead of the usual procedure of minimising the free 
energy is a well established fact of the replica method, appearing as soon as the number of replicas is 
less than 1 [2]. 

As a function of A , the thermodynamic potential has a maximum at: 

_ d(l - m) 1 _ 3 1 

where g* is the pair correlation of the liquid at the temperature T*. A study of the potential 

^{m, Ajnax), that equals (p{m), as a function of m then allows to find all the thermodynamic properties 
that we seek, using the formulas of the previous section. This step and the results will be explained 
below in sect. 8, where we shall compare the results to those of other approximations. 



7.3.3 Higher orders 

The systematic expansion of the thermodynamic potential i/j in powers of A can be carried out easily 
to higher orders. However the result involves some more detailed properties of the liquid at the 
effective temperature T* . For instance at second order one needs to know not only the free energy 
and pair-correlation of the liquid at temperature T*, but also the three points correlation. The results 
for the second order, that will be discussed in next section, will be obtained in the framework of the 
hypernetted chain molecular approach that is described in appendix. 
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7.3.4 Harmonic resummation 

One can obtain a partial resummation of the small cage expansion described above by integrating 
exactly over the relative vibration modes of the molecules. We shall use such a procedure here, that 
is a kind of harmonic expansion in the solid phase ^. 

We work directly with l/a = and start from the replicated partition function (201), within the 
quadratic expansion of the interaction potential v in the relative coordinates u^. (Clearly it is assumed 
that the l/a — > 0"*" limit has been taken, and that its effect is to build up molecular bound states). 
The exact integration over the Gaussian relative variables gives: 



i<j 

where the matrix M, of dimension 3N x 3N, is given by: 



Zm = J n d'^i exp -Pm ^ vin - r,) - ^^^TV log (/3M) (207) 



= Sij J2 '"i^'^^'^i ~ "^k) - v^iuin - rj) (208) 
k 

and v^y{r) = d'^v / dr^j^dry. We have thus found an effective Hamiltonian for the centers of masses rj of 
the molecules, that basically looks like the original problem at the effective temperature T* = T/m, 
complicated by the contribution of vibration modes that give the 'Trace Log' term. We expect that 
this should be a rather good approximation for the glass phase. Unfortunately, even within this 
approximation, it is not possible to compute the partition function exactly. The density of eigenstates 
of the matrix M is a rather complicated object and we have developed a simple approximation scheme 
in order to estimate it. 

We thus proceed by using a quenched approximation, i.e. neglecting the feedback of vibration modes 
onto the centers of masses. This approximation becomes exact close to the Kauzman temperature 
where m — ^ 1. The free energy is then: 

-^\og{m) log(27r) - J- log Z{T*) + (Tr log (/3M))* (209) 

Zm 2m m,N Zm 

that involves again the free energy and correlations of the liquid at the temperature T* . Computing the 
spectrum of M is an interesting problem of random matrix theory, in a subtle case where the matrix 
elements are correlated. Some efforts have been devoted to such computations in the liquid phase 
where the cigcnmodes are called instantaneous normal modes [60] . Here we use a simple resummation 
scheme that should be reasonable at high densities-low temperatures: it is described in the appendix 
. Using these results we can compute the replicated free energy Fm only from the knowledge of the 
free energy and the pair correlation of the liquid at the effective temperature T*. The results will be 
discussed in section 8. 



some loose sense the first order in tlie A expansion may be compared to Einstein approximation for the specific 
heat and the harmonic approximation we describe he may be compared to Debye's jellium. 

^^If the harmonic approximation were fully consistent, all the eigenvalues of H (the so called INN, Instantaneous 

Normal Modes [60]) should be positive. This is not the case, however the number of negative eigenvalues becomes very 
small at low temperature, still in the liquid phase, signaling that valleys can be approximately defined in this region. 
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7.3.5 Without replicas 

Up to now we have used the rephca theory. Rephca theory is a very powerful tool, but it has the 
disadvantage that many of the underlying physical hypothesis cannot be seen in a clear way. We will 
rederive some of the previous formulae without using the replica formalism [40, 43]. We only suppose 
that at low temperatures the phase space of the system can be approximately divided into valleys that 
are separated by high barriers [82, 83, 84]. In a first approximation each valley can be associated to 
one inherent structure [6, 85], i.e. one minimum of the potential energy and consequently there is a 
one to one correspondence among the valley at two different temperatures. 

Let us consider a system with N particles with Hamiltonian H{C), C denoting the generic con- 
figuration of the system. If we use the approach of the previous sections, the crucial step is the 
computation of the generalized partition function: 

Z{r,P) = exp(-iV7G(7;/3)) = ^exp(-7iV/„(/3)). (210) 

a 

Using the definition of the free energy in a valley 

/ dCexp(-7i7(C)) = exp(-7V7/(7,C)), (211) 

we obtain 

Z{r,P) = I dCexp(-7i7(C)-iV7/(/?,C) + 7V7/(7,C)) = 

I dC exp (-7if (C) - iV7/(/3, C) + iV7/(7, C)) , (212) 

where f{f3, C) = /(/?, C) — /(oo, C) and /(/3, C) is a function that is constant in each valley and it is 
equal to the free energy density of the valley to which the configuration C belongs. 

Before entering into the computation of /(/?, C) it is useful to make the so called quenched approx- 
imation, i.e. to make the following approximation inside the previous integral: 

exp(-A/(^, C)) = exp(-A(/(^))^), (213) 

where {fiP))^ is the expectation value of /(/3, C) taken with the probability distribution proportional 
to exp(— 7iJ(C)). The quenched approximation is exact if the temperature dependance of the energy 
of all the valleys is the same, apart from an overall shift at zero temperature. In other words we 
assume that the minima of the free energy have different values of the free energy but similar shapes. 
The quenched approximation would be certainly bad if we were using the free energy /(/3, C) at the 
place of f{(3,C) because the zero temperature energy strongly varies when we change the minimum. 
We finally find 

G{r, P) = Fl{i) + W) - /^(7), (214) 

being the free energy of the liquid {Sl{'^) and S^{j) are respectively the entropy of the liquid 
and of a valley). A simple algebra shows that 

S(7; 13) = Sl{i) - S,{^) + f'^{^) - f'^{f3), (215) 

where /' (/3) = df^{P)/d^. 
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In the liquid phase, we find out that the configurational entropy is given by 

= f3) = SUP) - Sp{(3). (216) 

The entropy of the hquid is the entropy of the typical valley plus the configurational entropy. 
The thermodynamic transition is characterized by the condition 

= 0. (217) 

In the glassy phase the free energy can be found by first computing the value of 7(/3) such that 

S(7(/3);/3)=0. (218) 

The quantity 7(/3) is the inverse of the effective temperature of the valley. It is easy to show (following 
[40]) that the previous formulae are completely equivalent to the replica approach. 

A strong simplification happens if wc assume that the entropy of the valley can be evaluated in 
the harmonic approximation where we only keep the vibrational contributions. For a system with M 
degrees of freedom the harmonic entropy of the valley near to a configuration C is given by 

S{m) = fin (^) - (ln(H(C7))) , (219) 

where 'H{C) is an M x M Hessian matrix (e.g. if H depends on the coordinates Xj we have that 
'^i,k = d'^H/dxidxk). The final result is just the same obtained with the replica method in the 
harmonic approximation if we put m = 7//?. 

8 The replica approach to structural glasses: some results 

8.1 Three approximation schemes 

We have seen up to now three approximation schemes. 

• The small cage expansion has been carried out directly to first order in section 7.3.2, and agree 
with the first order expansion within the molecular HNC approach. 

• The second approximation scheme is the harmonic resummation method. Again we have an 
explicit form (241) for the free energy per particle ^(m) only from the knowledge of the free 
energy and the pair correlation of the liquid at T* . Having this m dependance the procedure to 
get the thermodynamic results is entirely the same as that of the first order result. 

• The third approximation scheme is obtained by the expansion of the molecular HNC free energy 
to second order in the cage size, as described in the appendix. 

For each of the three approximation schemes mentioned above, we need to compute the free energy 
and the pair correlation of the liquid in a temperature range close to the glass transition. Wc will first 
consider three dimensional soft spheres [81] interacting through a potential v(r) = 1/r^^. We work 
for instance at unit density, since the only relevant parameter is the combination T = pT~^^'^. In this 
case e have used the hypernetted chain approximation to get both the correlation function g{r) and 
the free energy. Later on we will present the results for an LJ binary mixture within the harmonic 
approximation. 
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Figure 17: The effective temperature of the molecular liquid at the transition, T* = T/m* (left) and 
the quantity A/T, versus the temperature T, computed in an expansion to first order (dashed-dotted 
Une) and second order (full line) in the cage size A, and in the harmonic resummation (dashed line). 

8.2 Critical temperature and effective temperature 

We plot in fig. (17) the inverse of the effective temperature T*, equal to m*/T, versus the temperature 
T of the thermostat. The transition temperature is given by T* = T. This gives the ideal glass 
transition temperature. Within the first order expansion we find Tk — .14; the harmonic resummation 
gives Tk — .19 and the second order perturbation theory is Tk — .18 We see that the two best methods, 
the second order and harmonic resummation, are in good agreement and both give a critical value 
of r around T ~ 1.52. This value of T is in good agreement with numerical estimations of the glass 
transition of the soft sphere system, that range around 1.6 [86]. We also notice that the effective 
temperature stays relatively constant when the actual temperature varies. The effective temperature 
T* (that can be experimentally observed) is always near to Tk, independently from the value of the 
temperature T, 

8.3 Cage size 

In replica space the cage size characterizes the size of the molecular bound state, in the approximation 
of quadratic fluctuations, as defined in (194). Its physical meaning is easily established: In the glass 
phase at low temperatures one can approximate the movement of each atom as some vibrations in a 
harmonic potential in the neighborhood of a local minimum of the energy. The typical square size of 
the displacement is given by: 

A = {{n-{n)f) (220) 

that is the physical definition of the square size. The cage size is plotted versus temperature in fig. 17. 
The cage size is nearly linear in temperature, as it would be in a T-independent quadratic confining 
potential. This indicates that the local confining potential has little dependance on the temperature 
in the whole low temperature phase. 
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Figure 18: The free energy (left) and the specific heat (right) versus the temperature, computed in 
an expansion to first order (dashed-dotted line) and second order (full line) in the cage size A, and in 
the harmonic resummation (dashed line). In the right panel the dotted line is the specific heat of the 
liquid. 



8.4 Free energy, specific heat and configurational entropy 

If we plot the free energy versus the temperature one would see a strong consistency between the second 
order term of the small cage expansion and the harmonic resummation . Both data extrapolates at 
zero temperature to a ground state energy of order 1.95. This is related to the typical energy of the 
amorphous packings of soft spheres. More precisely, if we consider all the amorphous packings of soft 
spheres at unit density, we can count them through the zero temperature configurational entropy. The 
lowest energy where one can find an exponentially large number of such packings is the ground state 
energy of the glassy phase that we find equal to 1.95. However we have not taken into account the 
existence of a crystal: therefore we must first remove all crystal like configurations, i.e. configurations 
that correspond to a crystal with some local defects. These configurations can be characterized by the 
presence of delta functions at the appropriate values of the momenta. This procedure of identifying 
crystal like solutions has been explicitly done numerically in [82]. 

In fig. 18 (right panel) we plot the specific heat versus temperature. It is basically constant and 
equal to 3/2. This is the Dulong- Petit law (we have not included the kinetic energy of the particles, 
that would give an extra contribution of 3/2). This result is very welcome: if we had treated the 
crystal at the same level of approximation as we considered here for the glass, we would have the 
Einstein model for which the specific heat is also given by the Dulong-Petit law. Thus the specific 
heat of the glass is very near to that of the crystal, has it happens experimentally. Notice that it 
was not obvious at all a priori that we would be able to get such a result from our computations. 
The fact of finding the Dulong-Petit law is an indication that our whole scheme of computation gives 
reasonable results for a solid phase. 

In fig. 19 we show the configurational entropy versus the free energy at various temperatures, 
including the zero temperature case. We have included here for simplicity only the result from the 
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Figure 19: The configiirational entropy versus the free energy, computed within the harmonic 
resummation, at temperatures T = 0., .05, .1 (from left to right). 

harmonic resummation procedure. 

We notice that the various curves corresponding to different temperatures are not far from being 
just shifted one from another by adding a constant to the free energy. This indicates that the main 

effect of temperature is to add a constant (~ 3/2kT) in the energies of all amorphous packings. This 
correspond to the case where the vibration spectrum is approximatively state independent. 

8.5 The missing dynamical transition 

We Icnow that at the mean field level there exists a dynamical transition at a temperature Tjj larger 
than the thermodynamic transition temperature Tjc- This phase is characterized by the dynamic 
statement that a system will remain forever in the same valley, and its free energy is greater than 
the equilibrium one because it misses the contribution of the configurational entropy. This dynamic 
phase is just a mean field concept, that should disappear when corrections, such as activated processes, 
due to the short range nature of the potential, are taken into account. However if the barriers are 
sufficiently high, metastable states have a very large life time and they strongly affect the dynamics. 

In the framework of the harmonic resummation one finds that the approximation breaks down 
at small but positive e if the matrix of second derivatives has negative eigenvalues . From this 
point of view the appearance of negative eigenvalues signal the dynamic transition. Unfortunately in 
the chain approximation all the eigenvalues are positive at all temperatures and no dynamic phase 

^^'In the small cage expansion the perturbative method assumes that there is always a bound state: the breakdown of 
this assumption hardly be seen in a perturbative approach. 
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Figure 20: Left panel: analytical entropy of the liquid (upper line) compared with the numerical one, 
and analytical harmonic entropy (lower line) compared with the numerical results. The horizontal axis 
is r~'^. Right panel: analytical value of the complexity (upper line) compared with the numerical one 
(+ points), as functions of 



transition can be seen: the free energy is always well defined for small e.: the chain approximation 
may be reasonable at low temperature but it is certainly not good at high temperatures. One should 
use a better method to compute the spectrum, giving reasonable results also at higher temperature, 
e.g. following the approach of [68] (work in this direction is in progress). 

It is clear that a study of the dynamical phase transition should be done refining the tools than 
the one we have developed here. This is not surprising: the dynamical phase transition is present at a 
temperature higher than the static one and the approximations that we had used are especially good 
at low temperature. 

8.6 Lenhard-Jones binary mixtures 

There other model that is interesting to consider because they do not crystallize [87, 43]. Here we 
report the results fo is a realistic model for glasses, is given by a binary mixture of particles (80% large 
particles, 20 % smaller particles) interacting via a Lennard-Jones potential, introduced by Kob and 
Andersen [75]. This Hamiltonian should mimick the behaviour of some metallic glasses and it is one 
of the best studied and simplest Hamiltonian that do not lead to crystalization at low temperature. 
We first present the data at the density p = 1.2. 

The fact that the systems does not crystallizes implies that we can easily obtain information doing 
numerical simulations without serious difficulties [43]. Therefore one can use both an analytic or a 
numerical method to get information in the liquid phase. Let us start by presenting the numerical 
results [43]. A system oi N = 260 particles, in a cubic box with periodic boundary conditions at density 
p = 1.2 has been studied via Monte Carlo simulations. The entropy is obtained using the formula 
S{(3) = S{0) + d/3' {E{(3) — E{(3')). Given an equilibrium configuration the nearest minimum of the 
potential is found by steepest descendent and the computation of the 780 eigenvalues of 'H{C) does 
not present any particular difficulty. 

The results for the total entropy of the liquid and for the harmonic part are shown in fig. (8.6), 
left as function of T~'^ (a more detailed description of the simulations can be found in rcf. [43]). The 
entropy of the liquid is remarkable linear when plotted versus T^-^, as it happens in many cases [88]. 
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Figure 21: The specific heath coming from the x dependent part of the Hamiltonian as a function of 
the temperature 

In fig. (8.6), right, we show the configurational entropy as function of T~'^. We fit it with a 
polynomial of second degree in T~-^. The extrapolated configurational entropy becomes zero at a 
temperature Tk = -31 ± .04, where the error contains systematic effects due to the extrapolation 
(similar conclusions have been reached in ref. [83]). 

There are many methods to compute analytically the free energy in the liquid phase that lead to 
integral equations for the correlation functions. Unfortunately the simple HNC dos not works well. 
Here I present the results [43] obtained by mixing the HNC and MSA (mean spherical approximation) 
closures [89].. This technique allow us to compute the internal energy in the liquid with a reasonable 
approximation. The computation of the spectrum is more involved: it can be done with the same 
approach we used before that is described in the appendix. 

We can now put everything together [43]: the final analytic predictions for the liquid and harmonic 
entropy are shown in fig. (8.6), left. The predictions for the liquid entropy turn out to be very 
good, while there is a minor discrepancy for the harmonic entropy, likely due to the rather strong 
approximations in the analytic evaluation of the spectrum. The analytic configurational entropy 
is shown in fig. (8.6), left) It becomes zero at Tk = .32, that is our analytic prediction for the 
thermodynamic transition. 

We have in our hands all the tools to compute analytically the free energy in the low temperature 
case. If one computes the specific heat, one finds that the Dulong Petit law is extremely well satisfied 
in the low temperature region The value of 7(/?) weakly depends on /?: its value in the limit /? — oo 
is only about 10% higher that its value (i.e. /3x) at the transition temperature. 

Summarizing we have found a simple method that is able to use liquid theory method in the 
glasses phase putting in practice the old adage a glass is a frozen liquid. We are able to compute with 
a reasonable approximation the thermodynamics and with a little more effort we can compute the 
static and the dynamic structure functions. 

We may wonder what happens if we change the density p. The Kauzmann temperature should 

^^In the harmonic approximation the Dulong Petit law would be exact if we neglect the 7 dependence of Sj{f3). 
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be depend on p. A numerical estimated of the Kauzmann temperature can be done by studying the 
behaviour of the diffusivity: this leads to a numerical of a temperature Tq that would be near to T^. 
The results are shown in fig. 9 together with the analytic estimates for the same model [91] that are 
in reasonable agreement with the numerical data. 

9 Discussion and perspectives 

Within the equilibrium framework, we have implemented so far our general strategy using rather crude 
methods. These methods should be improved, and one should perform a more careful study of the 
molecular liquid. The extension to the case of hard spheres seem to be particularly interesting. 
There are also many other different problems that should be investigated, 

• A careful numerical study of the phase diagrams for coupled replicas would be welcomed in order 
to test the correctness of the detailed theoretical predictions 

• It would be very interesting to demonstrate numerically that the Kauzman transition exists 
beyond mean field theory and that the behaviour of the system is similar to the one we expect 
theoretically. Unfortunately known numerical techniques are not sufficient for thermalize a glass 
forming liquid at low temperature in a reasonable amount of CPU time: during the simulation 
the system is trapped in one of the exponentially large number of metastable states {0{exp{AN)) 
where ^ is a quantity of order 1),. The introduction of an appropriate lattice model, enough 
simple to be simulated in an very effective way, could be very useful in this respect. 

• Analytic and numerical method should be developed in order to compute the height of the 
activation barriers that should dominate the dynamics at low temperature (below the mode 
coupling transition). 

• The analytic computation of the spectrum of the instantaneous normal modes could be strongly 
improved. This could br used to study the properties of the Boson peak and to find analytically 
the value of the mode coupling and Kauzman temperature within the same approach. 

• The analytic approach to the computation of the thermodynamic properties should be extended 
to the case of hard spheres. 

• Least, but not the last, quantum glasses are a very wide territory that should be explore theo- 
retically with much more details. 

The study of the properties at equilibrium study is to be considered as a first step before dealing the 
dynamics (equilibrium and off-equilibrium). A very interesting and open problem is the computation 

of the time dependent correlation functions (and as a by-product the viscosity) in the region above Tk 
and below Tc- However a better understanding of activated processes in this framework is a crucial 
prerequisite. 

It is quite possible that in the next years we shall see progresses in some of the previously mentioned 
fields. 
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Appendices 

The probability distribution of the weights 

We need to compute the probability distribution of the weights w that are equal to: 

- 2Et^y , (221, 

where 

Z = 5^exp(-Ffe) . (222) 

k 

The probability of finding an energy Fj^ in the interval [F, F + dF] , is given by 

p{E) = exp(/?mF) , (223) 

and the possible values of the Fk go from — oo to +oo. 

This computation can be done with a brute force method [92] . One assumes that the index k goes 
from 1 to M and the F's belong to the interval [— oo, Fm]- One finally sends M to infinity, keeping the 
quantity Mexp(mFM) equal to a constant (the value of the constant is irrelevant). Esplicite formulae 
can be written and the appropriate approximations can be done in the limit M going to infinity using 
the saddle point limit. The computation is not too long, but it is not too simple. 

A much more clever method can be found in [93]. Here one tries to compute 



^ E ^k^ ■ (224) 



A;=l,oo 

without introducing the cutoffs M and Fm- After a few passages one finds that 



^^'^ = Jt7 . ■ (225) 

r(s)r(i - m) ^ ' 

In this way we can go backward and we finally find that the probability of finding a Wk in the interval 
[w,w = dw] is given by v{w)dw, where v{w) is given by eq. (39). 

The mean field approach for a ferromagnet 

It is interesting to see the techniques of Gaussian integration and saddle point at work in a case where 
we already know the result. The simplest case where this can be done is the infinite range ferromagnet 
[15], a model that can be exactly solved.. The Hamiltonian is given by: 

H=i C7{i)a{k)-h Y: ^(0, (226) 

i,k=l,N i=l,N 

where the sum over i, k is done over all the A'^(A'^ — l)/2 pairs of spins. An direct computation shows 
that in the limit N ^ oo the mean field approximation become exact and the magnetisation satisfies 
the mean field equation 

m = th(/3( Jm + h)) . (227) 
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We would like to prove this result using the saddle point method. 
At this end we can rewrite the Hamiltonian as 

^ = ^f E -h E (228) 

\i=l,Ar / /=l,Ar 



where we have neglected and addictive constant equal to — J, that originates from the sums of the 
terms with i = k. 

The argument runs as follows. In the limit N ^ oo, neglecting multiplicative constants that give 
no contribution to the free energy density) we can write that 



J2 f dmexp I + {13 Jm + ^h) ^ a{i) | = (229) 



dmexp{-Nf{m)) , (230) 
where 

/(m) = (3Jm^/2 - In (2ch(/3Jm + (3h)) . (231) 

Up to now everything was exact. The long range nature of the interaction and its homogeneity 
allows us to reduce the sum over variable to an one dimensional integral. This integral can be 
estimated by various means. In the limit N ^ oo we can us the method of the point of maximum. It 
gives 

Z ^ exp(-iV/(m*)) , (232) 
where m* is the minimum of the function f{m) and it is a solution of the equation 

1^ = . (233) 
am 

The previous equation coincide with the usual mean field equation eq. (227) . The correction to this 
result are proportional to 1/N (if we are not at the critical temperature) and they can be exactly 
computed by evaluating the corrections to the saddle point method. 

The spectrum of the instantaneous normal modes 

Here I would like to explain a fast method [94] for obtaining some analytic estimates on the spectrum 
of the matrix the matrix M, of dimension 3A^ x 3A^, given by: 

^{in)Uu) = Sij J2 ^M^^C^i - '^k) - v^u{ri - Tj) (234) 

k 

and Vni,{r) = d^v/dr^dri,. The fact that the diagonal terms of M do fluctuate complicate the analysis. 
This difficulty may be removed, if we notice that in this high density regime^^. there are many neigh- 
bours to each point, and thus a good approximation is to neglect the fluctuations of these diagonal 

■^*Hcrc and in what follows, wc have not written explicitly the density. We choose to work with density unity in order 
to simplify the formulae; however we assume that in the glsissy phase the density is sufficiently high, compare to the 
range of the forces, that an expansion in inverse powers of the density gives the correct result. At low density te approach 
of [95] gives the correct results. 
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terms and substitute them by their average value. We thus write: 

E - ^k) d'rg*{r)Av{r) = ro (235) 

k 

In this approximation (235) the diagonal matrix elements are all equal and can be factorized, 
leading to: 

(Tr log (M))* = m log( ro) + (Tr log (^5^,- ^ " ^^^i^-i^i " ^fc)) )* (236) 

This form lends itself to a perturbative expansion in powers of 1/ro that we assume to be a small 
number. The computation of the p-th order term in this expansion, 

= ( 5^ '^/^i/i2(''n ~ ^i2)-"'^Mp-iMp(^V-i ~ ''ip)^A'p/"i (''ip ~ ''n) (237) 

P'o 



II ...Ip 



still involves the p-th order correlation functions of the liquid at T*. We have approximated this 

correlation by introducing a simple chain approximation involving only the pair correlation. This chain 
approximation consists in replacing, for p > 2, the full correlation by a product of pair correlations. A 
most drastic approximation consist in selecting only those contributions that survive in the high density 
limit; systematic corrections can be computed following [94, 68]. Within the chain approximation, Tp 
is approximated by: 

Tp = j dxi...dxp g*{xi,....,Xp)[v^^^^{xi-X2).:V^^_^^^{xp-i-Xp)v^^^^{xp-xi)\ 

E y dxi...dxp [g*{xi - X2)v^^^j,^{xi - X2)\...[g* {xp - xi)v^^^j,^{xp - xi)] . (238) 



Ml— /Up " 



In this last form we need to compute a convolution that can be factorised through the introduction 
of the Fourier transform of the pair correlation function. We thus introduce the Fourier transformed 
functions a and b that are defined from the pair correlation g*{r) by: 

j d\ g*{r)v^,{T)e'^' = 5^, a{k) + - ^5^,) b{k) . (239) 

In terms of these Fourier transforms, the p-th order term in the 1/ro expansion is simply 

Tp = j d^k l^a{k) + '^b{k)J + (2) J d^k l^a{k) - ^b{k)J , (240) 

and the summation of the series over p is easily done and we finally find: 

(Trlog (/3M) =)*dlog(/3ro) 

+/d3, (aik)+m] ^^^^ (aik)-m\\ 



\2 



-I [ d'rg{r)Y^'-^ (241) 



where the function L3 is defined as: 

L3(x) = log(l-x)-Fx-Fa;V2 . (242) 
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The hypernetted chain approximation 

We derive here the form of the HNC free energy (252) for our molecular replicated system. One could 
use the standard diagrammatic method [96] , but here we shall follow the 'cavity' like method of Percus 
[97]. 

In the phase where the replica symmetry is broken, replicas arc correlated: it is convenient to 
consider N molecules with coordinates Xi,i G {1,..,N}. Each Xj stands for the coordinates of all 
atoms in molecule i: Xi = {x^}, a G {1, m}. The energy of the system is given by 

^ = E ^(^i' "^j) + E ^^^i) (243) 

i<j i 

where y(x, y) = X)a^(^" ~ v"")^ ^ intermolecular potential and the external potential f/(x) = 

uix"") has been introduced for future use. 

We shall need the following definitions. The one molecule density is 

^W = E(n'^«-^'^))' (244) 

i a. 

where the average (•) is done with respect to the Boltzmann measure exp(— /3iJ). 
The two molecules correlation {g) is defined as: 

p(2)(x,y) =^(n5«-x'')n5(4-^')) =p(x)<7(x,y)p(y) (245) 

i^j a b 

where we have also defined the pair correlation function g{-x.,-y), that goes to one at large (center of 
mass) distance. The connected pair correlation is: 

/i(x,y)=5(x,y)-l . (246) 

Functional differentiation gives: 

d^lpuly)) = " + P^'"^''^'"' ^^"^^^ ^^^^^ 

One can also introduce the direct correlation function c(x, y) through: 

d{-pui^)) _ 1 



dp{y) p(x) 



5(x-y)-c(x,y) . (248) 



The direct correlation is related to the connected pair correlation through the Ornstein-Zernike equa- 
tion c = (1 + hp)~^h that reads more explicitly: 

c(x,y) = /i(x,y) + j dxi/i(x,xi)p(xi)/i(xi,y) 

+ j (iXi(iX2^(x,Xi)p(xi)^(xi,X2)/9(X2)/i(X2,2/) + ... (249) 

The idea of Percus is to compute the pair correlation by considering the one point density with 
a molecule fixed at one point. Let us consider a problem where we have added one extra molecule, 
fixed at a point z = {z^ , ...,z^}. This extra molecule creates an external potential i7(x) = F(x, z). 
The one point density in the presence of this external potential, /^[/(x), is related to the density /9(x) 
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and pair correlation g(x, z) in the absence of an external potential through the conditional probability 
equation: 

/9c7(x) = p(x)5(x,z) . (250) 

The previous equations can be used to find a perturbative expansion of the logarithm of the 
correlation function. Doing the appropriate computations one find: 

log 5(x, z) + /3y (x, z) = Jdy c(x, y)/9(y)/i(y, z) . (251) 

Together with the inversion relation (249), this defines a closed set of equations for the one and two 
point molecular densities that are the HNC closure. In the real word the HNC equations are not exact: 
there are corrections to the r.h.s. 

In a similar way one finds that the free energy in the HNC approximation is a functional of the 
molecular density p(x) and the two point correlation ^(x, y). The result is: 

/3V' = ^y"c^X(iyp(x)p(y) [5(x,y)log5((x,y) -5((x,y) + 1 + /3i;(x,y)5((x,y)] 

- -^Tr f log(l + hp) -hp+ Ihphp) + - f dxp(x) log ^ (252) 

2m V 2 J m J e 

In the trace term all products are convolutions For instance the lowest order term in the small p 
expansion of the trace is: 

~\I ^^^y y)p(y)^(y^ z)p(z)^(z, x)p(x) (253) 

We would like to optimize the thermodynamic potential ip with respect to the molecular density 
p(x) and the two point function 5r(x, y). We shall work at low temperatures for which p should be 
nearly Gaussian. We thus choose an Ansatz for p of the type: 



(^)"'p--*'"'-e.p(-^i:(."-.V) (254) 



ab 



where the molecular density is parametrized by the single parameter A 
The ideal gas contribution (last term in (252) gives: 

J Yl d^x^pix) log ^ = N Q(l - m) log(27ry4) + ^(1 - m) - ^ logm - 1^ (255) 

The interaction term is more complicated, and it can be evaluated in the small cage regime. At 
the end of the day one can compute the correlation function g in the limit of small cage radius A, 
expanding in powers of A. In this one we recover the zeroth and the first order of the A expansion 
that we have obtained by a direct method. In the same way after a long computation [81] we obtain 
the second order in A that has been used in the main text. 



^^The first correction is proportional to and the proof of tliis statement can be done most easily using a diagram- 
matical approach; it is remarkable that a first order computation in h may be so accurate. 
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The trace could also be written as Tr(L3(— /ip), where L3 is defined in eq. (242). 
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